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We construct a Ginsburg Landau (GL ) theory to study the phases of liquid, sohd, superfluid, 
especially a possible supersolid and phase transitions among these phases in a unified framework. In 
this GL, we put the two competing orders between the solid component and the superfluid component 
on the same footing. We only introduce two parameters: i; which is the repulsive interaction between 
the normal component and the local superfluid mode and g which is a periodically changing chemical 
potential for the local superfluid mode. The microscopic origins of v and g are given. By using this 
GL action, we study the superfluid to supersolid transition, normal solid to the supersolid transition 
and analyze the conditions for the existence of the supersolid. The non-classical rotational inertial 
(NCRI) in the SS state is calculated and found to be isotropic in bcc and fee lattice, but weakly 
anisotropic in hep lattice. We study the elementary low energy excitation inside a supersolid. 
We find that there are one upper branch and one lower branch longitudinal "supersolidon" modes 
inside the SS, while the transverse modes in the SS stay the same as those inside the NS. We also 
determine the corresponding spectral weights of the two branches. We work out the experimental 
signatures of these elementary excitations in Debye- Waller factor, density-density correlation, vortex 
loop interaction and specific heat. The estimated excess entropy due to vacancies seems consistent 
with data measured in the specific experiment in Helium 4. Detecting the two supersolidon modes by 
various equilibrium and thermodynamic experiments such as X-ray scattering, neutron scattering, 
acoustic wave attenuation and heat capacity can prove or disapprove the existence of a supersolid in 
Helium 4. A toy model of supersolid wavefunction is analyzed. The difference and similarities with 
lattice supersolid are clearly demonstrated. Elementary excitations inside a Fulde - Ferrell-Larkin- 
Ovchinnikov (FFLO) state of superfiuid are also discussed. 



I. INTRODUCTION 



A solid can not flow. It breaks a continuous transla- 
tional symmetry into a discrete lattice translational sym- 
metry. There are low energy lattice phonon excitations in 
the solid. While a superfluid can flow even through nar- 
rowest channels without any resistance. It breaks a global 
C/(l) phase rotational symmetry and has the off-diagonal 
long range order (ODLRO) ^ . There are low energy super- 
fluid phonon excitations in the superfluid. A supersolid 
is a state which breaks both the continuous translational 
symmetry and the global U{1) symmetry, therefore has 
both the crystalline order and the ODLRO. The possibil- 
ity of a supersolid phase in *i/e was theoretically spec- 
ulated in 1970^"^. Andreev and Lifshitz proposed the 
Bose-Einstein condensation of vacancies as the mecha- 
nism of the formation of supersolid^. Chester wrote down 
a wavefunction which has both ODLRO and crystalline 
order and also speculated that a supersolid cannot ex- 
ist without vacancies or interstitialS.. Leggett proposed 
that solid '^He might display Non-Classical Rotational 
Inertial (NCRI) which is a low temperature reduction 
in the rotational moment of inertia due to the super- 
fluid component of solid '^He'^"^ . Leggett also suggested 
that quantum tunneling of He atoms between neighbor- 
ing sites in a crystal can also lead to a supersolid even 
in the absence of vacancies. Over the last 35 years, a 
number of experiments have been designed to search for 
the supersolid state without success. However, recently, 
by using torsional oscillator measurement, a PSU group 



lead by Chan observed a marked 1 ~ 2% NCRI of solid 
'^He at ~ 0.2K, both when embedded in Vycor glass 
and in bulk ^ife^. The authors suggested that the su- 
persolid state of '^He maybe responsible for the NCRI. 
The PSU group also detected about 10~* reduction in 
the rotational inertial in solid H^^, but by performing 
blocked cell experiments, the authors concluded that the 
reduction is classical due to the motion of HD impurities 
clustering. So the PSU group did not find any NCRI in 
solid hydrogen^. 

The PSU experiments inspired extensive both 
theoreticali^rJ^ and experimental^rii interests in the 
very intriguing supersolid phase of '^He. Now it was 
widely believed that disorders played key roles in the 
PSU experiments. There are two kinds of complemen- 
tary theoretical approaches. The first is the microscopic 
numerical simulatioiJ^. The second is the phenomeno- 
logical approacbi^^— . In this paper, we assume no disor- 
ders in our GL theory and study all the possible phases 
in different parameter regimes of the GL model. A solid 
has the order in the number density, while a superfluid 
has the order in the phase. The two phases are in to- 
tally different extremes of the state of matter. How to 
combine the two opposite extremes into a new state of 
matter, supersolid, is an important and interesting topic 
on its own. It is widely believed that the only chance to 
get a supersolid is that the solid is not perfect, namely, 
it is an incommensurate solid which has defects such as 
vacancies or interstitials, so the total number of bosons 
can fluctuate to give some rooms for the invasion of su- 
perfluid. Due to large zero point motions, there are in- 
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deed rapid exchanging between local '^He atoms in bulk 
'^He, but this local process will not lead to a global phase 
ordering. It was estimated from both X-ray scattering 
experiments^^ and some microscopic calculations^^ that 
the thermal excitation energies of a vacancy and inter- 
stitial are e„ ~ 10 K, a ~ 40 K, so thermal fluctuations 
favor vacancies over interstitials. Because both have very 
high energy, so the thermal generated vacancies and in- 
terstitials are irrelevant around 45 ~ 200 mK. However, 
it is still possible that there are quantum fluctuation gen- 
erated vacancies and interstitials even at T = 0. It is still 
not known if quantum fluctuations favor vacancies over 
interstitials or notl^. 

In this paper, we will construct a phenomenological 
Ginsburg Landau ( GL ) theory to study all the pos- 
sible phases and phase transitions in helium 4 system. 
We identify order parameters, symmetry and symmetry 
breaking patterns in all the phases. Particularly, we will 
address the following two questions: (1) What is the con- 
dition for the existence of the SS state ? (2) If the SS ex- 
ists, what are the properties of the supersolid to be tested 
by possible new experiments. Let's start by reviewing all 
the known phases in '^He. The density of a normal solid 
(NS) is defined as n{x) = no + n^e^^'^ = no + Sn{x) 
where n*^ = n_g and G is any non-zero reciprocal lattice 
vector. In a normal liquid (NL), if the static liquid struc- 
ture factor S(k) has its first maximum peak at fc„, then 
near km S{k) ~ ^ ^c(k'^-k'^ )^ • liquid-solid transi- 

tion is weakly first order, it is known that the classical 
free energy to describe the NL-NS transition is2i: 

+ Un ^ "'Gi"'G2"'G3"'G4^Gi-l-G2+G3-(-G4,0 +(-'^) 
Gi ,G2,G^^Gi 

where = r„ + c(G'^ — fep)^ is the tuning parameter 
controlled by the pressure or temperature. Note that 
the average density no does not enter in the free en- 
ergy. Because the instability happens at finite wavevec- 
tor, Eqn|T] is an expansion in terms of small parameter 
alone, it is not a gradient expansion ! The GL pa- 
rameters w and Un may be determined by fitting the 
theoretical predictions with experimental data. It is easy 
to see that EqnlT] is invariant under x — >■ x -f- a, n{x) — > 

n(x-l-a), n{G) — )• n(G)e*'^'° where a is any vector. In the 
NL, {n{G)) = 0, the translational symmetry is respected. 
In the NS, {n{G)) ^ 0, the symmetry is broken down 
to translational invariance under only a lattice constant 
a ^ R,G ■ R ^ 27m, n{x) n{x + R),n{G) -J> n(G). 
As shown in the Fig.l, the NL to NS transition only 
happens at finite temperature, so the classical theory is 
valid. Note that due to the lack of the Z2 symmetry of 
dn{x) — > —Sn{x), namely, ~^g^ there is always a 

cubic w term which makes the the NL to NS a 1st order 
transition. The w„ term which is invariant under the Z2 



symmetry is needed for the stability reason. 

Of course, the Superfluid to Normal Liquid transition 
at finite temperature in the Fig.l is the 3d XY transition 
described by^: 

U = K\V^\^ + t\tlj\^+u\ij\^ + --- (2) 

where ^ is the complex order parameter and t is the 
tuning parameter controlled by the temperature or pres- 
sure. Eqnl2]is invariant under the global U{1) symmetry 
4> "06*^ • In the NL, (-0) = 0, the symmetry is re- 
spected. In the SF, ^ 0, the symmetry is broken. 
The GL parameters K, t, u may be determined by fitting 
the theoretical predictions with experimental data. In 
this paper, we always assume u > 0. 

The coupling between n{x) and ip{^) consistent with 
all the symmetry can be written down as: 

/„t = gSn{x)\^j{x)\-' + v{5n{x)f\^{x)\^ + ■■■ (3) 

where Sn{x) — n{x) — uq — YI'g^^g^^'^^ ■ ^^te that the 
average density does not enter in the interaction ei- 
ther. Eqn|3]is an expansion in terms of two small param- 
eters Sn(x) and ipix). The • • • include the other terms 
in higher odd and even powers of Sn(x) which are sub- 
leading to the g and v term. In an effective GL theory, 
n{x) and ip{x) emerge as two independent order param- 
eters. 

Eqnl3] is invariant under both the translational sym- 
metry x — > X + a,n{x) — > n{x + d),ijj{x) — > ^{x + a) 
and the global U{1) symmetry — > -06*^- Note that it 
is important to keep both g and v term in the Eqn|3J 
because the g term changes the sign, while the v term 
is invariant under the Particle-Hole ( PH ) transforma- 
tion dn{x) — > —dn{x), so the sign of g makes a difference 
! Due to the two competing orders between the density 
fluctuation represented by {6n{x))'^ and the local super- 
fluid mode represented by we expect v to be always 
positive and is an increasing function of the pressure p. 
The positive v term is also needed for the stability rea- 
son. On the other hand, we can view 5 as a periodic 
chemical potential with average zero acting on -0. It is 
easy to see the coupling is attractive gy < for vacancies, 
but repulsive gi > for interstitials. From Eqn|3l we can 
classify three kinds of solids: If g = 0, the C-NS has 
the P-H symmetry, let's call this kind of PH symmetric 
C-NS as NS-PH. In general, g ^ 0, so there is no particle- 
hole symmetry in the G-NS, there are still two kinds: ( 
1) vacancy like C-NS where the excitation energy of a 
vacancy is lower than that of an interstitial, named NS- 
v. (2) interstitial like C-NS where the excitation energy 
of an interstitial is lower than that of a vacancy, named 
NS-i. We expect that in contrast to v, g is an intrinsic pa- 
rameter of solid Helium 4 which depends on the mass of 
a ^He atom and the potential between the ^He atoms, 
but not sensitive to the pressure p. In SS-v and SS-i, 
'0(x) stands for vacancies and interstitials respectively. 
The total density of the system is nt{x) = n{x) ± |V'(x)p 
for SS-i and SS-v respectively where n{x) is the density 
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of number of solid sites and |-!/'(a:)p is the superfluid den- 
sity. Both were treated as independent order parameters 
in the GL theory. 

The GL equations I1I2I3I are invariant under both the 
translational symmetry x — >■ a? + a,n{x) — > n{x + 

d),n{G) — > n{G)e^^'^,il;{x) — >■ ijjix + a) and the global 
U{1) symmetry — > ipe^^ . In a NL, — 0, = 0. 
In a SF, = 0, (V-) ^ 0. In a NS, ^ 0, (V') = 0. 
While in a supersolid, 7^ 0, (i/j) ^ 0. From the normal 
liquid ( NL ) side, one can approach both the solid and 
the superfluid. Inside the NL, t > 0, ip has a gap, so can 
be integrated out from EqnlSl we recover the solid-liquid 
transition tuned by in Eqnll]( Fig. 2 ). Inside the NL 
{n{x)) = no, the density fluctuations of Sn{x) is massive, 
so can be integrated out from Eqn|31 then we recover the 
NL to SF transition tuned by t in Eqnl2] ( Fig.l ). 

Quantum fluctuations can be incorporated by (1) 
n{x) — > n{x,T) and including ^pn{drn)'^ in EqnlT] 
(2) ip{x) — >■ %P(x,t) and including {x,T)drip{x,T) in 
Eqnl^l (3) Due to the lack of particle-hole symme- 
try in the normal solid, including additional terms like 
Sn{x,T)'ip^ {x,T)dri'{x,T) in Eqn. |3l We will explicitly 
consider the quantum fluctuations in section III where 
we will discuss the quantum phase transition from the 
superfluid to the C-NS ( either NS-v or NS-i ). How- 
ever, as shown in this paper, the quantum fluctuation 
terms are very important in the zero temperature quan- 
tum phase transition from SS-v ( SS-i ) to NS-v ( NS-i ) 
driven by the pressure ( see Fig. 2 ) and to determine the 
low energy excitation spectra in a given phase. 

The rest of the paper is organized as the following: 
In sec. II, for the first time, we construct a quantum 
Ginsburg-Landau action to study SF to the NS transition 
and also explicitly establish the Feynman relation from 
the QGL. In Sec. Ill, taking into account the couplings 
between the n and ip sector encoded in Eqnl3l we study 
the superfluid to supersolid transition which is a simul- 
taneous combination of the SFDW transition in ip sector 
driven by the roton condensation at fco — kr and the NS 
transition in the n sector driven by the divergence of the 
structure function at ko ~ kr = kn discussed in Sec. II. 
We also sketch the global phase diagrams to be confirmed 
and analyzed in the following sections. In Sec. IV, we ap- 
proach the SS phase from the NS side and discuss the 
SFDW in SS-v and SS-i respectively. We analyze care- 
fully the conditions for the existence of SS-v. We explic- 
itly show that the SS-v is the possible ground state when 
the gv is sufficiently negative ( Fig. 3). The SS-i is less 
likely, but we still analyze SS-i in the same footing as 
SS-v. We also classify several common SS-i lattice struc- 
tures. In Sec. V, by renormalization group analysis, we 
study the universality class of zero temperature quantum 
phase transition from normal solid (NS ) to supersolid 
(SS) driven by the pressure. In Sec. VI, we calculate 
the NCRI's in both the SF and SS states and find the 
NCRI in the hep SS lattice may be weakly anisotropic. 
In Sec. VII, we work out the non-topological elementary 
excitations and corresponding spectral weights inside the 



supersolids in both the isotropic solid case and the hep 
lattice structure case. Then in Sec. VIII, we study topo- 
logical elementary excitations, namely, vortex loops and 
vortex lines in the SS by performing a duality transforma- 
tion to the vortex loop representation. We also estimate 
the very large vortex core size and low critical velocity 
in the SS state. Then in the following sections, we study 
the experimental signatures of these low energy elemen- 
tary excitations: In Sec. IX, we make key predictions on 
the elastic X-ray scattering amplitudes from all SS-v and 
SS-i structures classified in Sec. IV by calculating the 
Debye- Waller factor and the density-density correlation 
function. In Sec. X, we study the specific heat in the 
SS. In Sec. XI, inspired by the results achieved in Sec. 
VII, we discuss the elementary excitations inside a Fulde 
- Ferrell-Larkin-Ovchinnikov (FFLO) state of supercon- 
ductors. Conclusions are summarized in the final section 
XII. In the appendix A, we discuss the properties of a 
tight-binding toy SS-v wavefunction. In the appendix B, 
we compare the properties of the SS with those of SS on 
extended boson Hubbard model in a lattice. Two short 
reports of these results appeared ini^ii^. 



II. SUPERFLUID TO NORMAL SOLID 
TRANSITION 

Let us start with the SF to the NS transition. In the 
superfluid state, if the multi-quasiparticle part can be 
neglected in the dynamic structure factor, the Feynman 
relation between the Landau quasi-particle dispersion re- 
lation in the ^ sector and the equal time structure factor 
in the n sector holds: 

^(9) - IT^, (4) 

In the (7 — >■ limit, S{q) ^ q, ijj{q) ^ q recovers the 
superfluid phonon spectrum near g = in Fig. lb. The 
first maximum peak in S{q) corresponds to the roton 
minimum in uj{q) in the ^^2 sector ( Fig. lb ), namely, 
kn = kr- As one increases the pressure p, the interaction 
u in Eqn|3] also gets bigger and bigger, the first maxi- 
mum peak of S{q) increases, the roton minimum gets 
smaller and smaller^^. Across the critical pressure^ = pc, 
there are the two possibilities. (1) The resulting solid is a 
commensurate solid, then (ip) = 0. In the NS, there is no 
remanent of the roton inside the SF, the supersolid phase 
does not exist as an equilibrium ground state. This is the 
SF to the C-NS transition. As to be shown in section IV, 
this happens when \g\ is sufficiently small ( Fig.l ). (2) 
The resulting solid at high pressure is an incommensurate 
solid with zero point quantum fluctuations generated va- 
cancies or interstitials whose condensation leads to the 
formation of the SS-v and SS-i respectively. In section 
IV, this happens when \g\ is sufficiently large ( Fig. 2 ). 
We will discuss case (1) in this section, then the most 
interesting case (2) in the next section. 
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The effective action inside the SF is : 

C[5n,e] = iSndrO + ]^Ps{VOf + ]^5nVn{q)5r 



(5) 



where ps is the superfluid density and Vn{ci) = a — bq^ + 
cq^ with a > 0, 6 > is the density-density interaction 
between the ^ife atoms. 

In the SF state, it is convenient to integrate out 5n in 
favor of the phase field 9 to get the phase representation 

where one can identify the compressibihty — 
\\mVn{q)\q^Q = a. The dispersion relation of the super- 
fluid modes including higher orders of momentum can be 
extract ed^i: 



= [PsVnmo' - '^Psq^ia - bq^ + cq^) 



(7) 



where one can see the superfluid phonon velocity — 

Ps/n- 

It is easy to see that the dispersion relation indeed 
has the form shown in Fig. lb with a roton minimum. 
Because the original instability comes from the density- 
density interaction Vn{q) , it is convenient to integrate 
out the phase field in favor of the density fluctuation 
operator 6n. Neglecting the vortex excitations in 9 and 
integrating out the 9 in EqnjS] leads to: 

C[Sn] - hni-q, -|- Vn{q)]Sniq, cj^) (8) 

where we can identify the dynamic pseudo-spin 
density-density correlation function Sniq,i-^n) = 

2 

{Sn{-q,~ujn)Sn{ci,ujn)) = ^2ff(^)^2 where 
= psVniq) is the " momentum dependent " 



spin wave velocity. Obviously v (q = 0) 



Ps/k. 



From the pole of the dynamic density-density cor- 
relation function, we can identify the speed of sound 
wave which is exactly the same as the spin wave ve- 
locity. This should be expected. From the analyti- 
cal continuation iLo„ uj + id and taking the imagi- 
nary part, we can identify the dynamic structure factor: 
S>{q,uj) = S„{q)S{LL> - v{q)q) where S„{q) = psq7T/2v(q) 
is the equal time density correlation function shown in 
Fig. lb. As q — > 0, Sn{q) — > q- The sum rules are: 



du>S'^{q,u>) = psqTr/2v{q) = Sniq) 



dujujS^{q,uj) = psq^TT/2 

) 

dujS>{q,uj)/uj = ps7T/2v^{q) 



(9) 



where the first one just defines the equal time density 
correlation function, the second one gives the / sum rule 
which can be used to define the superfluid density ps, 



the third one can be used to define the compressibility k: 
dujS>{q,u})/uj\g^o = k-k/2. 
The Feynman relation which relates the dispersion 
relation to the equal time structure factor is: 



uj{q) 



J^^ dujujS>{q,uj) ^ p.nq'^ 



(10) 



which takes exactly the same form as Eqn|4] if we iden- 
tify ps ~ 1/m with m the mass of ^He mass. Therefore, 
we recovered the Feynman relation from our GL theory 
Eqnl5]which gives us the confidence that Eqnl5]is the cor- 
rect starting action to study the the SF to SS transition. 
The density representation EqnjS] is dual to the phase 
representation EqnlHl However, the phase representation 
Eqn|6] contains explicitly the superfluid order parameter 
■0 ~ e*^ which can be used to characterize the superfluid 
order in the SF phase. While in Eqn|8l because the or- 
der parameter tp is integrated out, the superfluid order 
is hidden, the signature of the superfluid phonon mode 
is encoded in the density sound mode, so it is not as ef- 
fective as the phase representation in describing the SF 
state. However, as shown in the following, when describ- 
ing the transition from the ES to the NS, the density 
representation EqnIS] has a big advantage over the phase 
representation. 

Because the instability is happening at g = go instead 
of at q — 0, the vortex excitations in 9 remain uncritical 
through the SF to SS transition. Integrating them out 
will generate interactions among the density 6n: 



C[Sn] = ^Sn{-q, 



-a;„)[ — ^ + Vn{q)]Sn{q,Un) 



PsQ^ 



w{Sn)^ + u{Sn) 



(11) 



where the momentum and frequency conservation in the 
quartic term is assumed. Note that the (w„/q)^ term 
in the flrst term stands for the quantum fluctuations of 
6n which is absent in the classical NL to NS transition 
EqnlT] Because of the lack of Z2 exchange symmetry, 
there is a cubic term in Eqn llll Expanding Vn{q) near 
the roton minimum qo leads to the quantum Ginsburg- 
Landau action to describe the SF to the NS transition: 



C[Sr 



1 



6n[Aniol + r + c{q^ - ql'f]dr 



w{Snf + u{5n)'^ 



(12) 



where 



Pel - p and Ap 



1 

Pa 95 



which is non-critical 



across the transition. Just like EqulU because the insta- 
bility happens at the finite wavevector g = go, EanlT2lis 
not a gradient expansion, but an expansion in terms of 
the small order parameter Sn. Again the average den- 
sity no does not appear in Eqn ll2l The generic transi- 
tion driven by the collapsing of roton minimum is from 
SF to NS instead of from the SF to the supersolid ( 
SS ). In the SF, r > 0, (V-) ^ 0, {6n) = 0, In the NS, 



r < 0, (0) = 0, {6n) = Ed 



JG-S 



G '-G'- 



where G are the 
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(a) 



FIG. 1: (a) The theoretical phase diagram of GL model 
EansfTl2land[3l in P versus T plane. This phase diagram only 
happens when \g\ is sufficiently small, so the normal solid is 
a C-NS ( Fig.3 ). T controls thermal fluctuations, while p 
tunes quantum fluctuations. SF is the superfluid phase, C- 
NS is the commensurate normal solid phase which still does 
not have the P-H symmetry, it is likely to be a vacancy-like 
C-NS ( NS-v ). NL is the normal liquid phase. The super- 
solid phase in Fig. 2 is absent. The {^pi) and ('02) stand for 
the condensations in zero momentum and the flrst recipro- 
cal lattice vectors respectively. Thick ( thin ) lines are 1st 
( 2nd ) order transitions. The critical temperatures Tsf of 
NL to SF transition drops slightly as the pressure p increases 
because of the quantum fluctuations. The SF to the C-NS 
transition is described by the QGL Eanll2l (b) The struc- 
ture factor and the separation of low ( phonon ) and high ( 
roton ) momenta regime in the dispersion relation in the SF. 
The Feynman relation between the two curves are explicitly 
derived in EgnllOl 



shortest reciprocal lattice vector of the resulting lattice. 
The corresponding phase diagram for the SF to the NS 
transition is shown in Fig.l. 



III. SUPERFLUID TO SUPERSOLID 
TRANSITION 

In the last section, we discuss the SF to the NS tran- 
sition with the order parameter Sn. In this section, we 
discuss the SF to the SS transition with both order pa- 
rameters ■0 and Sn. Then we have to also include the 
couplings between Sn and ■0 sector encoded in EquIS] 
Starting from the SF side with increasing pressure, we 
develop the theory based on the two facts (1) there is a 
roton minimum in the superfluid state (2) The instability 
to solid formation is driven by the gap diminishing at the 
roton minimum. As shown in this last section, the fact 
(2) is guaranteed by the Feynman relation Eqn|4] which 
relates the roton minimum to the peak of the structure 



factor. Across the phase boundary p = Pci in Fig. 2, the 
resulting solid could be an in-commensurate solid ( IC- 
NS )with vacancies or interstitials even at T = whose 
condensation leads to to the formation of the SS-v and 
SS-i respectively where (tp) ^ Q^^tl^. There is still some 
remanent of the SF in the IC-NS, the supersolid phase 
does exist as an equilibrium ground state. In this section, 
we assume \g\ is sufficiently large in Fig. 2 and study the 
SF to the SS transition across p = Pci- Because ■0 is 
also critical through the transition, we can not simply 
integrate it out like in the last section. In fact, in the 
effective GL theory, we have to treat n and ip on the 
same footing. From Eqn[T]and Eqnl2l we can see that n 
and Ip have very similar propagators, so the lattice for- 
mation in n sector with n{x) = no -I- YI'g '^g^*'^'^ where 
G = kn and the superfluid density wave ( SFDW ) forma- 
tion in -0 sector with {ip{x)) = -00 + e'^^ J2m=i ^mS^^'"'^ 
where (5,„ — kr — kn ^ li^ja ~ 27r/3.17A have to 
happen simultaneously. From Hansen- Verlet criterion^i^ 
when S{kn)/no is sufficiently large, solidification in the n 
sector occurs, so the roton minimum remains finite just 
before its condensation. The SFDW p = \tp\'^ is sim- 
ply locked to ( or commensurate with ) the underlying 
normal solid ( n ) lattice. In fact, this locking is dic- 
tated by the p density- Sn density couplings in EqnlHl 
If the coupling g is attractive < 0, the SFDW just 
coincides with the n lattice. If it is repulsive gi > 0, 
then it simply shifts the SFDW by suitable constants 
along the three unit vectors in the direct lattice. These 
constants will be determined in the next section for dif- 
ferent n lattices. Namely, the supersolid states consist 
of two inter-penetrating lattices formed by the n lattice 
and the ip2 superfluid density wave. In fact, in a care- 
fully prepared super-pressured sample, the roton mini- 
mum survives up to as high pressure as Pr ~ 120 bar. 
This fact suggests the roton minimum in the meta-stable 
state in a super-pressured sample is replaced by a SFDW 
which is commensurate with the n lattice in the stable 
SS state. Strikingly, this pressure ^ 120 bar is close to 
Pc2 ~ 170 bar in Fig. 2 where the NCRI was extrapolated 
to vanish in the PSU's experiments-. This consistency 
also lead some support to the above picture. Obviously, 
when the pressure is so high that p > Pc2 in Fig. 2, C-NS 
is the only possible ground state, any remanent of the SF 
completely disappears. So the vertical axis in the Fig. 2 
shows the SF-SS-NS series as the pressure P increases at 
T — Combining the roton condensation picture^'^ in the 
last section with the results reviewed in the introduction, 
we can sketch the following phase diagram Fig. 2 of the 
complete QGL Eqns. [H [l]and[l 



As can be seen from the Fig. 2, starting from the SF 
side, as the pressure is increased at a given temperature, 
there are two possible states (1) At very low T < Tss, 
the Bose-Einstcin condensation ( BEG) of tp2 leads to the 
supersolid state where {ipi) ^ 0, {tp2) ^ 0, (ng) ^ 0. The 
instability happens at a finite wavevector k = kr instead 
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T.P -p 

(a) 

FIG. 2: (a) The theoretical phase diagram of GL model 
Eans lll2l andl3lin the pressure P versus temperature T plane. 
This phase diagram only happens when \g\is sufficiently large 
( Fig. 3 ). The SS is likely to be the vacancies induced super- 
solid ( SS-v ). Tss is an effective measurement of the strength 
of the coupling constant (/ 7^ in Eqn|3l The {ipi) and (?/)2) 
stand for the condensations in zero momentum and the first 
reciprocal lattice vectors respectively ( See Eans ll7l and [TS] 
). Thick ( thin ) lines are 1st ( 2nd ) order transitions. The 
SF to the SS transition is a simultaneous combination of the 
SFDW transition in sector driven by the roton condensa- 
tion at fco = kr and the NS transition in the n sector driven 
by the divergence of the structure function at ko — kr — kn- 
(b) The structure factor and the separation of low ( phonon 
) and high ( roton ) momenta regime in the SF. 

of at fc = 0. (2) At higher temperature Tss <T < Tsf, 
there is a direct SF to NS transition. The SS state is cer- 
tainly different from a conventional normal solid phase 
where (V'l) = {1P2} = 0, (n^) 7^ 0. In addition to the 
conventional translational and rotational orders charac- 
terized by by n, the SS also has the ODLRO character- 
ized by -01 and When decreasing the temperature 
at a given pressure, Hp < Pc, the NL becomes SF at 
T = Tsf, the instability happens at fc = 0. li p > Pc, the 
NL becomes a NS first at T = T™, then there is a SFDW 
onset transition from the NS to a SS phase at T = T55. 

Because the SF to SS transition driven by the roton 
condensation can be either weakly or strong first order, 
in principle, Eqn ll2l works precisely only in the SF side, 
it is not easy to study the precise nature of the SS state 
from the SF side. It turns out that it is more convenient 
to study the properties of the SS state from the NS side 
in the next section. 



IV. THE NORMAL SOLID TO THE 
SUPERSOLID TRANSITION 

Starting from the NS side, we assume that the local 
tunneling processes in the NS-PH only leads to local 



fluctuations of ^p with a gap A(p) in Fig. 3, so there is 
no long-range phase coherence and no supersolid in this 
case ( Tss = ). Taking the NS-PH state as the refer- 
ence state, we will show that if the coupling g in Eqn|3]is 
sufficiently negative ( or positive ), the supersolids SS-v 
or SS-i can be realized by adding small number of va- 
cancies to NS-v or interstitials to the NS-i ( Fig. 3). The 
attractive ( repulsive ) interaction g is shown to be cru- 
cial to raise the normal solid to the vacancy ( intersti- 
tials ) induced supersolid transition temperature Tss-v 
( Tss-i ) above the zero temperature in the Fig. 3. In 
fact, the temperature Tss becomes an effective measure 
of the coupling strength g. We find that SS-v is more 
likely than SS-i. However, in order to be complete, we 
study both cases on the same footing. It is also construc- 
tive to compare SS-v with SS-i even though the SS-i is 
unlikely to be relevant to the Helium 4 system. When 
the two kinds of SS show different properties, we treat 
them differently, when they share the same properties, we 
treat them just in the same notation SS. The SS phase 
naturally and consistently fits into the phase diagram. 

In the NS, ip stands for the vacancies or inter- 
stitials coming from the large zero point quantum 
fiuctuations^!^!^^. Inside the NS, the translational sym- 
metry is already broken, we can simply set dn{x) — 

n{x) — uq = X]'g "'g'S**^^ put it into Eqn|3]to look 
at the effects of the coupling constants g and v. Imag- 
ine that at a given pressure p > Pci ~ 25 bar, if tuning 
5 — >■ 0, but keeping v intact, then the normal solid be- 
comes an asymptotically P-H symmetric normal solid ( 
NS-PH ) in Fig.3. Inside the NL, the mass of V' at fc = 
is t — T ~ TxY which sets up the temperature scale of 
the problem. Inside the NS-PH, it is easy to see that the 

mass of ^ at fc = is ijvs-PH = t + ''^YI'g 
we take the temperature scale t = T — Txy as the ref- 
erence scale, then t^s-PH = T + A.{p) where A(p) = 
«E'g - Txy > is the T = gap for the local 

superfluid mode at fc = in the NS-PH. Because as 
the pressure p increases, the repulsive interaction v also 
increases, so it is reasonable to assume that A(p) is a 
monotonic increasing function of p. Because it is a first 
order transition across Pci , just like the roton gap > 
remains finite just before the first order transition, A(p) 
also remains finite just after the first order transition, 
namely, A(p+ ) > 0( Fig.3 ). 

Taking this NS-PH as a reference state, we then grad- 
ually turn on g and see how the ground state evolves. In 
the presence of the periodic potential of n[x) lattice, V' 
will form a Bloch wave. In principle, a full energy band 
calculation incorporating the interaction u is necessary 
to get the energy bands of ^. Fortunately, qualitatively 
important physical picture of GL Eqnsll] [21 [3] can be 
achieved without such a detailed energy band calculation. 
We expect that the self interaction u only leads to some 
renormalization of the parameters t and K . Writing the 
g term in momentum space g'^^ n{G)ip* {G)'ip{0) where 
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quantum fluctuations may favor vacancies over intersti- 
tials. However, for simplicity, we assume fi{gi) is also a 
monotonic increasing function of gi , so we can discuss va- 
cancies and interstitials induced supersolids on the same 
footing. In the following, we discuss SS-v and SS-i re- 
spectively. 



Vacancies induced supersolid: SS-v 



FIG. 3: (a) The gap A(p) of the local fluctuating superfluid 
mode tp in the P-H symmetric commensurate normal solid 
( NS-PH ) which exists only at g = 0. We assume it is a 
monotonic increasing function of the pressure p. (b) The zero 
temperature phase diagram of g^ verse the pressure p. The 
NS-PH only exists at g^ = 0. Any gv 7^ will transfer the 
NS-PH into the NS-v. If the experimental is along the path 
I, then it is a direct 1st order SF to NS-v transition in Fig.l, 
if it is along the path II, it is a 1st order SF to SS-v, then a 
2nd order SS-v to NS-v transition in Fig. 2. 



G 7^ and integrating out ^{G), we get a perturbative 
expansion on the eigen-energy e^(0) at the origin K = Q, 
it is easy to see the n' s term is of the form —(—5)". Up 
to the third order gjj, the series is: 



£^(0) = tNs-PH - glP 



n{Gi)n{G2)n{-Gi - G2) 
K f^G\K ^G\ 



where the term linear in g vanishes, fi = v,i stands for 
vacancies and interstitials respectively and the shortest 
reciprocal lattice vector is G = 27r/a with a ^ 3.17A the 
lattice constant of the solid '^He, the n{K) ~ e^^'° 
where K is any reciprocal lattice vector can be taken as 
a Gaussian where a is the width of the Gaussian. 

For vacancies, < 0, without writing out the coeffi- 
cients explicitly, the expansion is 



e„(0) = t NS-PH 



9l 



\9v\ 



\9v\ 



(14) 



so the coefficient has the same sign. Assuming the series 
converges, for any g^ < 0, we can write e«(0) = t — fv{gy) 
where the fv{gv) > 0, /i,(0) = is a monotonic increasing 
function of g^ and likely has no upper bound. 

For interstitials, gi > 0, without writing out the coef- 
ficients explicitly, the expansion is 



ei(0) = tNS- 



PH 



(15) 



so the coefficient has oscillating sign. Assuming the se- 
ries converges, we can write 6^(0) — t — fi{gi) which holds 
for any gi. Because of the oscillating nature of the ex- 
pansion coefficients, it is hard to judge the nature of the 
function of fi{gi) except we know /i(0) = 0. The dif- 
ferent expansion series of fv{gv) and fi{gi) indicate that 



The mass of ipy was decreased to t^,^, = T + A(p) — 
fv{gv) = T - Tss-v where Tss-v{p) fv{gv) - M.p) ( 
Fig. 3). Because fv{gv) is a monotonic increasing function 
of l^^l and /«(0) = 0, defining a critical value fv{gvc) = 
A(pci), then when < l^^d, fv{gv) < A(pci), the 
■01, mode remains massive, namely (ipy) = 0. The C- 
NS remains to be the ground state even at T = 0. It 
is important to stress that even the solid is a C-NS, it 
still does not have the P-H symmetry. For g^ < 0, the 
C-NS is a vacancy-like C-NS ( named as NS-v ) where 
the vacancy excitation energy e„ is lower than that of the 
interstitial Ci ( Fig. 3b ). 

If \9v\ > \9vc\, then Tss-v{Pci) = Uigv) - A(pci) is 
raised above the zero temperature, the SS-v state ex- 
ists in the Fig. 3b. Tss~viPci) is also proportional to the 
superfluid density measured in the experiments. The re- 
sulting solid is an in-commensurate solid with vacancies 
even at T = whose condensation leads to (ipv) ^ 0. Of 
course, the number of vacancies is quite small. The 
SS-v state has a lower energy than the NS-v state at suf- 
ficiently low temperature. As the pressure increases to 
Pc2, Tss-v{Pc2) = fv{gv) - A(pc2) = ( Fig.3 ). Then 
fvigv) = A(pc2), so Tss-vip) = MPc2) - A(p) which be- 
comes an effective experimental measure of the energy 
gap A{p) in Fig. 3a. Indeed, if reflecting the A{p) in 
Fig. 3a with respect to the horizontal axis and then shift 
by A(pc2) = fvigv), then one recovers the Tss-vip) in 
Fig. 2. In the following, substituting the ansatz {ipi{x)) — 

ae^^i and {M^)) = e'^^ A^e''?™'-^ where Q„ - Q 

into Eqnlll we study the effects of n lattice on -0 = 
"01 + V'2- From Eqnl31 we can see n{x) acts as a periodic 
potential on ip. In order to get the lowest energy ground 
state, we must consider the following 4 conditions: (1) 
because any complex V' ( up to a global phase ) will lead 
to local supercurrents which is costy, so we can take ip to 
be real, so Qm have to be paired as anti-nodal points. P 
has to be even (2) as shown from the Fcynman relation 
EqnUl QrajiTi = 1, • • • ,P are simply P shortest recip- 
rocal lattice vectors, then translational symmetry of the 
lattice dictates that e{K = 0) = e{K = Qm), ipi and ■02 
have to condense at the same time (3) The point group 
symmetry of the lattice dictates A^ — A and is real (4) 
for the vacancies case, the interaction is attractive g^ < 0, 
from Eqnini the SFDW-v simply sits on top of the n lat- 
tice, so the Superfluid Density wave p = |'0P simple sits 
on top of the n lattice as much as possible. This is reason- 
able, because vacancies are hopping on near the lattice 
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FIG. 4: The phase diagram of T versus g at a given pressure 
Pel < p < Pc2- The finite temperature transitions denoted by 
the dashed hne II at a given > \gvc\ { or Qi > gic ) is in 
the 3D XY universahty class described by Ean l22l However, 
if \gv\ < \gvc\, then the ground state at T = is a vacancy-hke 
C-NS( named as NS-v ). Similar thing can be defined for an 
interstitial-like C-NS ( named as NS-i ) when gi < gic- 



sites. From EqnlSl the attractive interaction also favors 
4>{x = R/2) ~ where R/2 stands for any interstitial 
sites which are in the middle of lattice points. It turns 
out that the the 4 conditions can fix the relative phase 
and magnitude of V'l a-nd "02 to be 02 = ^i, A = a/-P, 
namely: 



2 

— ^ cos Qr. 

ni—l 



x) 



(16) 



where ^po = ae'^ depends on the temperature and pres- 
sure. Note that in contrast to a uniform superfluid, the 
magnitude of ip is changing in space. This field satisfies 
the Bloch theorem with the crystal momentum fc = 
and the Fourier components are ip(K = 0) = a,ip{K = 
Qm) = a/ P. They have the same sign and decay in mag- 
nitude. In principle, higher Fourier components may also 
exist, but they decay very rapidly, so can be neglected 
without affecting the physics qualitatively. 



B. Interstitial induced supersolid: SS-i 

The discussion is quite similar to the vacancy case ex- 
cept that (1) \i gi < gic, the C-NS is an interstitial-like 
C-NS ( named as NS-i ) where the interstitial excita- 
tion energy is lower than that of the vacancy Cy (2) if 
gi > gic, the resulting solid is an in-commensurate solid 
with interstitials even at T = whose condensation leads 
to (ipi) 7^ 0, the SS-i state exists in the Fig. 2. The ar- 
guments to determine the lattice structure of the SS-i 
goes the same as those in the vacancies case except the 
condition (4), for the SS-v, the interaction is attractive 
gv < 0, the SFDW-v simply sits on top of the n lattice. 
However, for the interstitials case, the interaction is re- 
pulsive gi > favors tp{x = 0) ~ 0, so the Superfiuid 
Density wave p = can avoid the n lattice as much 
as possible. This is reasonable, due to the competing 



of the the two orders, the superfluid component emerges 
from the places where the normal solid component is sup- 
pressed . It turns out that the the 4 conditions can fix 
the relative phase and magnitude of ipi and ip2 to be 
^2 = ^1 + TT, A = a/P, namely: 



^o(l 



2 

— cos Qr. 

m—1 



■x) 



(17) 



where ■i/'o = cie*^ depends on the temperature and pres- 
sure. Note that the crucial sign difference from the va- 
cancies case which make an important difference from 
the X-ray scattering from SS-v and SS-i to be discussed 
in the next section. Again in contrast to a uniform 
superfluid, the magnitude of -0 is changing in space. 
This field satisfies the Bloch theorem with the crys- 
tal momentum fc = and the Fourier components are 
ip{K = 0) = a,ilj{K = Qrn) = -a/P. They oscillate in 
sign and decay in magnitude. In principle, higher Fourier 
components may also exist, but they decay very rapidly, 
so can be neglected without affecting the physics quali- 
tatively. 

In the following, we will discuss different lattice struc- 
tures of n when P = 2, 4, 6, 8, 12 respectively. Obviously, 
the Superfluid Density wave p = has the same Bra- 
vais lattice structure as the the n lattice. However, as 
shown in the following, even rt is a Bravais lattice, the 
SFDW may be the same Bravais lattice plus a few basis. 

(a) P=2: Qi ~ — Q2 — Q are a pair of anti-nodal 
points. They are the two shortest reciprocal lattice vec- 
tors generating a 1 dimensional lattice embedded in a 3 
dimensional system. The fleld is ^'(2?) = a(l — cos Q ■ x). 
There is a superfluid density wave formation transition 
inside the normal solid which is a 2nd order transition 
in the universality class of 3D XY model( Fig. 2). The 
local SFDW operator p[ = \^(x)\'^ — a'^(l ~ cosQ ■ x)^. 
It breaks translational invariance only along 1-dimension 
which is similar to Smectic-A or Smectic-C phase in the 
liquid crystaP^. The maxima of the SFDW ipmax = 2a 
appear exactly in the middle of lattice points at a = ^Si. 
So the SFDW forms the dual lattice of the Id lattice 
which is also a Id lattice. 

^(b) P=4: Q3 = -Qi,Qi = -g2,Qi ■ Q2 ^ 0, 
Qi, 1 = 1,2,3,4 form the 4 corners of a square. They are 
the four shortest reciprocal lattice vectors generating a 2 
dimensional square lattice embedded in a 3 dimensional 
system. The field is '0(^') = a[l — ^{cosQi-x+cosQ2-x)], 
The local SFDW operator p[ = \il!{x)\'^. The maxima of 
the SFDW 'ipmax — 2a appear exactly in the middle of 
lattice points at a = i(ai -f 02). So the SFDW forms the 
dual lattice of the square lattice which is also a square 
lattice. 

(c) P=6: there are two cases need to be discussed sep- 
arately, (cl) Qi,i = 1,2,3,4,5,6 form the 6 corners of 
a hexagon. They consist of the 6 shortest reciprocal lat- 
tice vectors generating a 2 dimensional triangular lat- 
tice embedded in a 3 dimensional system. The field is 
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ip{x) = a[l — ^(cosQi ■ X + C0SQ2 • X + cosQs • x)]. The 
maxima of the SFDW i^max = 3/2a appear in the middle 
of lattice points at a = ±i(ai +02). They form the dual 
lattice of the triangular lattice which is a honeycomb lat- 
tice. The honeycomb lattice is not a Bravais lattice which 
has two triangular sublattices A and it can be consid- 
ered as one triangular lattice A plus a basis. This can be 
understood intuitively: there are two equivalent ways to 
shift the n lattice, one way to get the sublattice A, the 
other to get the sublattice B. Putting A and B together 
forms the SFDW which takes the honeycomb lattice (c2) 
Qi,i = 1,2,3,4,5,6 are the 6 shortest reciprocal lattice 
vectors generating a cubic lattice. The maxima of the 
SFDW ipmax = 2a appear exactly in the middle of lat- 
tice points at the 8 points a = ^{ztai ± 02 ± 03). So the 
SFDW forms the dual lattice of the cubic lattice which 
is also a cubic lattice. 

(d) P — 8: Qi,i — I,-- - ,8 form the 8 shortest re- 
ciprocal lattice vectors generating a bcc reciprocal lattice 
which corresponds to a fee direct lattice. The field is 
tp{x) = a[l — j{cosQi-x+cosQ2-x+cosQ3-x+cosQ4-x)]. 
The maxima of the SFDW ipmax = 2a appear in all the 
edge centers such as (1/2, 0, 0) etc. and the centers of any 
cube such as (1/2, 1/2, 1/2) etc. It is easy to see these 
points can be achieved by simply shifting the n lattice 
by (1/2, 1/2, 1/2), so the SFDW also forms a fee direct 
lattice. 

(e) P = 12: = 1, • • ■ , 12 form the 12 shortest 
reciprocal lattice vectors generating a fee reciprocal lat- 
tice which corresponds to a bcc direct lattice. The field 
is ip{x) = a[l — i(cosQi ■ x + cos (52 ■ x + cosQa ■ x + 
cos Qi ■ X -f cos Q5 ■ X -\- cos Qe ■ x)]. The maxima of the 
SFDW Ipmax = 4/3a appear along any square surround- 
ing the center of the cube such as (1/2, /3, 0) or (1/2,0,7) 
etc. In fact, one can achieve the SFDW lattice by shift- 
ing the the center of bcc to the 3 face centers, so all these 
points are on the edge centers and face centers which are 
only discrete points on the square surrounding the center 
of the cube. We expect the continuous whole square is 
due to the artifact of the approximation ip(x = 0) ^ 
imposed. So the SFDW ^He in Vycor glass takes a bcc 
lattice. 

Unfortunately, a spherical k — Q surface can not lead 
to lattices with different lengths of primitive reciprocal 
lattice vectors such as a hep lattice. This is similar to the 
classical liquid-solid transition described by EqnlT] where 
a single maximum peak in the static structure factor can 
not lead to a hep latticeSi. Another difficulty with the 
hep lattice is that the hep lattice is not a Bravais lattice, it 
consist of two inter-penetrating simple hexagonal lattices 
shifted by a = ^ai + ■^(12 + ■^a^. Here we can simply take 
the experimental fact that n forms a hep lattice without 
knowing how to produce such a lattice from a GL theory 
Eqnll] Despite the technical difficulty, because for an 
idea hep lattice c/a = an hep lattice has 12 nearest 

neighbors, so its local environment may resemble that 
of an fee lattice. We expect the physics ( except the 



anisotropy of NCRI in the hep lattice to be discussed in 
the next two sections ) is qualitatively the same as that 
in fee direct lattice. From the insights achieved from 
the other lattices, one can achieve the SFDW lattice by 
shifting the hep lattice of n by a = |(ai -I- 02) + -jaa. 



C. Discussions on both SS-v and SS-i 

In fact, for both SS-v and SS-i, we can write n and 
Tp sector in a more symmetric way: n{x) = uq + 

Y^gnge^^'^ ,ip{x) — ipo + YI'g'^g^^^^ ■ ^^^^ *° 

see that in the SS, the translational symmetry x — > 
X -\- a,n^ — >■ n^e^'-^'"' ,ip^ — >■ ip^e^'~^'°' is broken down to 

a = R,ng ^ n^e^^'^ = ng^ipg — >■ ipge^^'^ = ip^^. In 
the Fig. 2, in the NL side, as the temperature is lowered, 
the symmetry breaking happens in Tpi at k = 0, the NL 
gets into the SF. However, as shown in this section, in 
the NS side, the symmetry breaking happens in both the 
ipi &i k — Q and ip2 sector at P discrete points simultane- 
ously, the NS gets to the SS state at a much lower critical 
temperature Tss- The results achieved from the NS side 
in this section indeed confirm Fig. 2 achieved from the 
roton condensation picture in the SF phase in the last 
section. 



V. THE ZERO TEMPERATURE TRANSITION 
FROM NS TO SS DRIVEN BY THE PRESSURE 

The analysis so far focused on finite temperature and 
mean field level. Some interesting physics near the finite 
temperature NS to SS transition in Fig.l was explored 
in^'* by considering the coupling of elastic degree of free- 
doms to the SF mode. For example, the sound velocity 
will acquire a dip similar to the specific heat cusp in the 
A transition in superfludi Helium. In this and following 
sections, I will push the Ginsburg-Landau (GL) theory in 
the previous sections to zero temperature and to include 
all the possible low energy fluctuations above the mean 
filed solutions achieved in the previous sections. Partic- 
ularly, I will work out the new elementary low energy 
excitations including vortex loop excitations in a SS and 
study how they defer from the low energy excitations in 
solids and superfluids. In principle, if these elementary 
low energy excitations can be detected by X-ray scatter- 
ing, neutron scattering, acoustic wave attenuation and 
heat capacity experiments in solid Helium 4 can prove or 
disprove the existence of the supersolid in Helium 4. In 
practice, the detection may still be complicated by sam- 
ple quality. No matter if a supersolid indeed exists in 
Helium 4, these results should be interesting in its own 
and may have application in other systems. 

So far, we only look at the mean field solutions corre- 
sponding to vacancies and interstitials. Here we discuss 
excitations above the mean field solutions. It turns out 
that the excitations in both cases are the same, so we 
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discuss both cases at same time. In the SS-v and SS-i, 
the wavefunctions can be written as 

ipss ^'4'oil± -"^cosQn,- x), Vo = IV'o|e'^ (18) 

m— 1 

where ± sign corresponds to SS-v and SS-i respectively. 
Obviously, there are also topological defects in the phase 
winding of 9 which are vortices. At T ^ Tss, the vor- 
tices can only appear in tightly bound pairs. However, 
as T — Tgg, the vortices start to become liberated, this 
process renders the total NCRI to vanish above T > T55. 
In addition to the superfluid 9 mode in SS states, there 
are also lattice phonon modes u in both n sector and ip2 
sector. However, it is easy to see that the coupling EqnI2] 
is invariant under x — >■ x + u, n{G) — ^ ri(G)e**^ ", V'(G) 
4>{G)e^^'^ , so the lattice phonon modes in tp are locked 
to those in the conventional n lattice. This is expected 
because there is only one kind of translational symmetry 
breaking, therefore only one kind of lattice phonons. In- 
side the NS side, the translational symmetry is already 
broken, so we can parameterize the density deviation or- 
der parameter Sn{x, t) = n{x, r)— no and the SF complex 
order parameter Tp{x^T) as: 

/ 

5n{x,T) = ^nge*'5-(5'+«(5'^^)) 

G 

^(f,T) = ^o(^,T)[l±iX]e^'5-(-"+"(^~'^»] (19) 

G 

where the iPq{x,t) = |'(/'o(a;, r)|e'^^^''^-' is the SF order 
parameter, u{x, r) are the 3 lattice phonon modes, the 
± means vacancy or interstitials induced supersolids re- 
spectively, n*g — the " / " means the sum over the 

shortest non-zero reciprocal lattice vector G and P is the 
number of them. From Ean ll91 we can identify the SS 
density order parameter pq{x^t) = q^G-u{x,t)^ r^-^^ effec- 
tive action to describe the NS to SF transition at T = 
consistent with all the lattice symmetries and the global 
C/(l) symmetry is: 

+ O^li^Ualjlpldrlpo + a]^pUal3\'lpo\'^ ^ (20) 

where r — p — pc2 with pc2 ^ 170 bar ( Fig. 5), p„ is 
the normal density, — ^{daUp + dpUa) is the strain 
tensor, XaP'yS are the bare stress tensor dictated by the 
symmetry of the lattice, it has 5 (2) independent elastic 
constants for a hep ( isotropic ) lattice. For a uniaxial 
lattice such as hep lattice, all the coefficients Cq.^, a"^, a^^ 
take the same form Cap — Cznanip+c_\_{5ai3 — nanp) where 
n is a unit vector along the uni-axisi^i^. In the NS state 
r > 0, {tpo^XjT)) = 0, the 3 lattice phonon modes u{x,t) 
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FIG. 5: Possible temperature T versus pressure p phase di- 
agram of Helium 4. SF is the superfluid phase, SS is the 
supersolid phase, NS is the normal solid phase, NL is the nor- 
mal liquid phase. T^dxy is the 3d XY transition. Tm is the 
1st order melting transition. The dot is the zero temperature 
transition from the NS to the SS which is a transition with 
mean field exponents z = 2,1^ — 1/2, r; = 0. Thick (thin) fine 
is 1st (2nd) order transition. This figure can be obtained from 
Fig.2 by rotating 90°. 



become the 3 ordinary ones. While inside the SS state 
r < 0, {ijjo{x,T)) 7^ 0. From the parameterizations of 
'iJj{x,t), we can see if the prefactor (7/10(2?, t)) 7^ 0, then 
■^(x, r) condenses at both G — and any other non- 
zero reciprocal lattice vectors G to form the superfluid 
density wave ( SFDW ) = \ip{x,T)\'^ inside the SS. The 
a^^ and a;i^^ couplings come from the original couplings 
Sn{x,T)tp^ drip and Sn{x,T)\'tp\'^ respectively in the GL 
in Sect. I. If setting all the couplings between tJjq and Ua 
vanish, the ipo sector describes the SF to Mott insulator 
transition in a rigid underlying lattice^S. So we will 
study how the NS to SS transition at T = 0,p = Pc2 in 
a rigid lattice is affected by its coupling to a quantum 
fluctuating lattice. Under the Renormalization group ( 
RG ) transformation, r' = x' = x/b and ip' = ip/Z. 
If we choose z ^ 2,Z = b-"^/"^, the g' = gb'^''^. It is 
well known that the SF to Mott insulator transition in 
a rigid underlying lattice has the mean field exponents 
with z = 2, 1/ = 1/2, 77 = at d > 2^6. We also choose 
u'a — Ua/Z, then = b~^pn, so the lattice phonon 
kinetic energy term is irrelevant near the QCP. It is easy 
to see Oq = 6~''/2-iQ_g^ go ug is always irrelevant, a'l = 
so both g and ai's upper critical dimension is 
du = 2, so, in principle, a e = 2 — d expansion is possible 
for both g and oi. However, both are irrelevant at d = 3. 
We conclude that NS to SS transition at T = remains 
the same as that in a rigid lattice described by 

£t=o = ipldripQ + Capdaipldptpo + r\ipa\'^ + g\ipo\'^ -I 

(21) 

Namely, it is a transition with mean field exponents z = 
2,i. = l/2,r, = 0. 

If neglecting the r dependence by setting Ua{x,T) — 
Ua{x),ipo{x, t) = ipoix), then Eanl20l reduces to the clas- 
sical action studied i n^^i^^ . For the classical case, x' = 
x/b,i)' = ip/Z,u'^ = Ua/Z, if we choose Z = b^^''^^^, 
then g' = gb'^"'^, a[ = b^~'''/^ai, so both g and ai's upper 
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critical dimension is = 4. So in principle, a e = 4 — d 
expansion is possible for both g and Oi, the putting e — I 
for d = 3. In^'^, it was shown that due to the specific heat 
exponent of the 3d XY model a = —0.012 < 0, the a;!^,^ 
coupling is irrelevant, so the NS to SS transition remains 
to be a classical 3d XY transition at finite temperature 
given by: 

Ct = ifjvslVV^oP +ijvs|^o|' + UNs\H^ + ■■■ (22) 

where Ins = T - Tss and Tss-vip) = fvidv) - A(p) = 
A{p,2) - A(p) and Tss-^{p) = M9^) ~ A(p) = A{p,2) - 
A(p). Obviously, it is still a 3d XY transition. It was 
shown ini^ that the coupling to the phonon mode u will 
not change the universality class at finite temperature. 
As the pressure increased to Pc2, Tss-v or Tss-i are 
suppressed to zero, the system becomes a C-NS where 
("0) = 0. It is important to stress that even at p > pc2, the 
solid is a C-NS, it still does not have the P-H symmetry 
because g ^ 0. For g^ < 0, the C-NS is vacancy like 
( NS-v ) where the vacancy excitation energy is lower 
than that of the interstitial. For gi > 0, the C-NS is 
interstitial like ( NS-i ) where the interstitial excitation 
energy is lower than that of the vacancy. 



VI. THE NCRI OF SF AND SS STATES AT 
T = AND T > 0. 

In order to calculate the superfluid density ps explic- 
itly, we need to look at how the system's free energy 
responses to a fictitious gauge potential A. We find that 
when A » Ac, Ps{T = 0) = p = J d'^x\^{x,T)\^ = 
J d''a;(|V'i(a?, t)P -I- \ip2{x , t)\'^) where the crossing terms 
between and 1/12 drop out due to the momenta conser- 
vation. 

In the SF state, at low T, the quantum fiuctuations 
induced by the pressure are important. Let's first look 
at the quantum phase fluctuations. The phase fluctu- 
ation action is given by Cp = j^^T,^u:„ 1 10i{uJn + 

k^)\9{k, aj„)p where 9 — j- controls the strength of quan- 
tum phase fluctuations and the superfluid phonon veloc- 
ity has been set equal to 1 for simplicity. It is easy to 
see that at T = 0, {9'^{x,t))t=o is infra-red ( IR ) fi- 
nite, so the quantum phase fluctuations alone will not 
lead to any instability. However, it will lead to super- 
fluid density depletion even before reaching the phase 
boundary of SF to SS transition in Fig.l and Fig. 2, al- 
though the depletion may be quite small. This fact 
explains why Tsf{p) bends to the left slightly as the 
pressure p increases. At finite T, the thermal fiuctua- 
tions (6i2(f,T))T - (6'2(f,T))T=o - T-^'^ lead to ps{T) = 
Ps{T = 0) — cT^ at d = 3. It is well known the superfluid 
density p{T) ~ p(T = 0) - aT^, while the Bose con- 
densation density nh{T) nb{T — {)) - bT^. So strictly 
speaking, ip sector can only describe the Bose condensa- 
tion. This is expected, because the n sector in the SF 



phase also contributes to the superfluid density, but not 
to the BEC. 

Then let's look at the roton fluctuations. Setting the 
roton gap A2 = A^, at T = 0, the quantum roton 
fluctuations {(I)\{x,t))t=o ^ log -^^ where A is the up- 
per cutoff is IR logarithmic divergent as A which 
signifies the instability to the lattice formation. Due 
to this IR divergence, the 1st order SF to SS transi- 
tion may happen well before A becomes zero, namely, 
at A = Ac > 0. This is consistent with the picture 
described in Sec. III. At finite T, the thermal roton fiuc- 
tuations {4>2{x,t))t — {4>\{x,t))t=o ^ (logA)e~T when 
T< A< A. 

In the SS-v and SS-i states, the n{x) forms a lat- 
tice, at the same time, the unstable roton part is re- 
placed by a stable SFDW formation commensurate with 
the underlying n lattice. Obviously, the n{x) normal 
lattice takes away the vast majority of density from 
the superfluid density even at T = 0. It can be 
shown that superfiuid density from the "01 sector is 
isotropic pi ~ Ka?, while the superfiuid density from 
"02 sector turns out to anisotropic in hep lattice p2.ij ~ 

Z]m=i \^mVQmiQm]/Q'^ wherc A„ = ±a/P for SS-v or 
SS-i. Therefore, the total superfiuid density in SS phase 
is: 

1 P 

m— 1 

(23) 

Taking P = 6, the anisotropy is quite small. Solid ^He 
in a bulk takes a hep lattice with c/a ^ 1.63 which is quite 
close to the idea value c/a = ^8/3. The three primitive 
reciprocal lattice vectors are Gi = G2 = -^2-- , G3 = ^ . 
We can estimate the anisotropy of the NCRI in the 
hep lattice. If the rotation axis is along the c axis, 
the NCRI is pii - Ka" + VrYL=i \^^?QmiQ^ilQ'' 
If the rotation axis is along the a ( or 6 ) axis, the 

NCRI is P33 - Ka" + vr \Am\^Qra3Qm3/Q^ ■ The 

anisotropy mainly comes from the ip2 sector. Setting 
s = G1/G3, for the idea value s = > 1, so pn > P33. 
Namely, the NCRI response is larger when one is rotat- 
ing the sample around the c axis than that when one is 
rotating the sample around the a or 6 axis. However, as 
the pressure is increased, s decreases, the anisotropy of 
the NCRI also decreases. In the PSU experiments, the 
samples are poly-crystal, the relative orientation of the 
rotation axis to the c axis is not known, so it's hard to 
test this prediction with poly-crystals. 



VII. NON-TOPOLOGICAL ELEMENTARY 
EXCITATIONS AND SPECTRAL WEIGHTS IN 
THE SS 

Classical non-equilibrium hydrodynamics in SS was in- 
vestigated for a long time^i^. These hydrodynamics will 
break down at very low temperature where quantum fine- 
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tuations dominate. However, the quantum nature of the 
excitations in the SS has not been studied yet. Here, 
we wiU study the quantum characteristics of low energy 
excitations in the SS. For example, how the phonon spec- 
tra in the SS differ from that in a NS and how the SF 
mode in the SS differs from that in a SF. Inside the SS, 
{^o{x,t)) = a, we can write ^q{x,t) — ^Ja + Spe'^^'^^'^^ 
and plug it into the Egnl^ 

c = iSpdre + pi^dc^edpO + hpSo^Sp 

+ al^pUaiiidrO + a]^pUai}5p^ (24) 

where we already dropped idrd term which is irrelevant 
inside the SS phase ( although it is very important in 
describing the T = NS to SS transition discussed in 
the last section ), 5*0 is the bare density-density correla- 
tion function. Integrating out the massive magnitude 5p 
fluctuations, we get: 

S 

C = — ^{idrO + al^pUapf + p%pda9dp9 

+ alpUapiadrO + ■ ■ ■ (25) 

Expanding the square, we get the effective action describ- 
ing the low energy modes inside the SS phase: 

c = ^[K{dr0)^ + pipd^edpd] 

where k is the SF compressibility defined in EqnlS] and 
/o^^ is the SF stiffness which has the same symmetry as 

0°^, o-ap = a"^ -I- Soa\p where So{k,uj) is the bare SF 
density correlation function. Again, the p„ is the normal 
density, the XajS^s is the stress tensor. Obviously, the last 
term is the crucial Berry phase coupling term which cou- 
ples the lattice phonon modes to the SF mode. The factor 
of i is important in this coupling. By integration by parts, 
this term can also be written as aapidrUisdaO + dTUadpO) 
which has the clear physical meaning of the coupling be- 
tween the SF velocity daO and the velocity of the lat- 
tice vibration drUfj. It is this term which makes the low 
energy modes in the SS to have its own characteristics 
which could be detected by experiments. The invari- 
ance under the Galilean transformation^^iS dictates that 

Oq,^ — Pn5a,p — Pa,l3- 

In this section, we neglect the topological vortex loop 
excitations in EqnlMl In the next section, we will discuss 
these vortex loop excitations in detail. In the following, 
we discuss two extreme cases: isotropic solid and hep 
lattice separately. Usual samples are between the two 
extremes. 



A. Isotropic solid 

A truly isotropic solid can only be realized in a highly 
poly-crystalline sample. Usual samples are not com- 
pletely isotropic. However, we expect the simple physics 
brought about in an isotropic solid may also apply qual- 
itatively to other samples which is very poly-crystalline. 

For an isotropic solid, Xa/s-yS = XSa/sS-ys + p-iSa-yS/ss + 
SasS/3-y) where A and fi are Lame coefficients, p^ p — 
P'''Sa,p,aa.p = a6a,p where a ^ pn - Ps- fn w„) space, 
Eqnl24| becomes: 

^is = i[p„a;^ + (A + 2/i)q2]|u,(g',a;„)p 

+ \[^i^l + Psq^Mq,^n)\^ 

+ \[pn^l + pq^]\M<i,^n)\^ (27) 

where ui{q,uJn) = iqiUi{q,ujn/ q is the longitudinal com- 
ponent, ut{q,u!n) = i^ijqiUj{q,uJn) I q are transverse com- 
ponents of the displacement field. Note that Egnim 
shows that only longitudinal component couples to the 
superfluid 9 mode, while the two transverse components 
are unaffected by the superfluid mode. This is expected, 
because the superfluid mode is a longitudinal density 
mode itself which does not couple to the transverse 
modes. 

From Eqn l271 we can identify the longitudinal- 
longitudinal phonon correlation function: 

, X ^ i^^l + Psq^ 

^""'"'^ {k^I + Psq^){pnwl + (A + 2p)q^) + aWu^l 

(28) 

The {69) and {uiO) correlation functions can be simi- 
larly written down. By doing the analytical continuation 
iujn — > uj + i5, we can identify the two poles of all the 
correlation functions at Wj. = v'^q^ where the two ve- 
locities v± is given by v\ = [k{X -f 2/i) -f psPn -I- ± 

(k(A -I- 2p) + psPn + aP'Y - AnpsPniX + 2p)]/2npn. If 
setting a = 0, then v\ reduces to the longitudinal phonon 
velocity vf^ = (A -I- 2p)/ pn and the superfluid velocity 
v1 = Ps/ K, respectively. Of course, the transverse phonon 
velocity vf^ = /J,/p„ is untouched. For notation simplic- 
ity, in the following, we just use Vp for vip. Inside the SS, 
due to the very small superfluid density ps, it is expected 
that Vp > Vg. In fact, in isotropic solid iJe^, it was mea- 
sured that vip ^ 450 — 500to/s, vt ^ 230 ^ 320to/s and 
Vs ~ 366to/s near the melting curve^^. It is easy to show 
that v+ > Vp > Vs > w- and v\ + v'^ > v^ + v^, but 
v+V- = VpVs, so v+ + V- > Vp + Vs ( see Fig.l ). Note 
that because the Galilean invariance dictates a — p„ — p^, 
for Ps <^ Pni one can see PsPn + a? ^ PsPn, so vj^ {v- ) 
are considerably above ( below ) Vp { Vs ), so the two su- 
persolidons, especially the softening of the lower branch, 
may be easily distinguished by possible neutron scatter- 
ing experiments. 
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argued in^ that there must be a diffusion mode of va- 
cancies in the NS. Inside the SS, the vacancies condense 
and lead to the extra superfluid mode. So the diffusion 
mode in the NS is replaced by the SF mode in the SS. The 
analysis in this subsection is similar to the atom-photon 
polaritons of atoms inside a cavity2Zi2^ or exciton-photon 
polaritons of excitons in quasi-two dimensional quantum 
wells inside a planar cavity^S. 



FIG. 6: (a) The elementary low energy excitations inside a 
supersolid. The coupling between the phonon mode ujp = Vpq 
( the upper dashed line ) and the superfluid mode lUs = Vi,q { 
the lower dashed line ) leads to the two new "supersolidon " 
modes ij± = v±q ( solid lines ) in the SS. Their corresponding 
spectral weights are listed in Ean l29l and 1301 These two new 
supersolid modes and the spectral weights should be detected 
by in-elastic neutron scattertings. (b) The excitations inside 
a lattice CDW supersolid discussed in the appendix B. It has 
one gapped mode with a CDW gap and a gapless superfluid 
mode. The corresponding spectral weights are worked out 
in^. 



By doing the analytical continuation iw„ — > cj + zJ, we 
can take the imaginary part and find: 

-[d{uj- v+q) - d{u} + v+q)\ 



- vt 2pnV+ q 



1 



w+ - vi 2pnV- q 



[6{uj - v^q) ^ 6{uj + v^q)] (29) 



It is easy to see that the second term can be achieved 
from the first term just by i;+ -f-^- w_ . Setting a — 0, then 
v+ = Vp, V- = Vs, the second term just vanishes, the 
first term recovers the excitation spectrum of the lattice 
phonons. When a 0, then Eqnl29] becomes a mixing 
of the lattice phonons and superfluid phonons, the first 
and second term give the excitation energies and the two 
corresponding spectral weights. 
Very similarly, we can find 



[S{uj - v+q) - 6{uj + v+q)] 
-[S{uj - v-q) - 5{uj + v-q)] (30) 
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It is easy to see that the second term can be achieved 
from the first term just by o u_ . Setting a = 0, then 
f+ = fp, f_ — Vg, the first term just vanishes, the second 
term recovers the excitation spectrum of the superfluid 
phonons. When a ^ 0, then Eqn l30l become a mixing of 
the lattice phonons and superfluid phonons, the first and 
second term give give the excitation energies and the two 
corresponding spectral weights. 

The two longitudinal modes in the SS can be under- 
stood from an intuitive picture: inside the NS, it was 



B. hep crystal 

Usual single hep crystal samples may also contain dis- 
locations, grain boundaries. Here we ignore these line 
and plane defects and assume that there are only vacan- 
cies whose condensation leads to the superfluid density 
wave inside the supersolid discussed in Sect. IV. 

For a uni-axial crystal such as an hep lattice, the action 

is: 

^^hcp = ]^[pn{drUaf + Kii{ul^+uly)+2Ki2U^xUyy 

+ KSSUI^ + 2Kl3{UxX + Uyy)Uzz 

+ 2{Kn-Ki2)uly + Ki4ul, + ul,)] 



+ [ttzdzU^ + a±{dxUx + dyUy)]idr 



(31) 



If q is along z direction, namely q^ 0,qx = qy = 0, 
then Eqn l31l simplifies to: 

^hcp = ^[Pn^l+ K33ql]\Uz{qz,UJn)\'^ 

+ l[i^^l + plql¥{qz,^n}\^ 

- iazqzU>nUz{-qz,-(^n)0{qz,UJn) 

+ ^[p„0jl+KMql/4:]\ut{qz,i0n)\^ (32) 

where \ut{qz,uJn)\'^ = |ua:('7z, t^n)!^ + \uy{qz,ujn)\'^ stand 
for the two transverse modes with the velocity = 
Ki^i/Apn- The superfluid mode only couples to the longi- 
tudinal Uz mode, while the two transverse modes Ux,Uy 
are decoupled. EqnI321is identical to EqnI27]after the re- 
placement Uz — >■ It;, ifas — A -I- 2/x, Oz a. It was found 
that vip ^ bAOm/sjVt ~ 250m/ s when q is along the z 
direction-^^. Fig. 2 follows after these replacements. 

Similarly, we can work out the action in the xy plane 
where — 0,qx 7^ 0, qy ^ 0. Then Uz mode is decoupled, 
only u,j.^Uy modes are coupled to the superfluid mode: 



^^cp = l^[Pn{drUaf + Kii{ul^ + uly) 

XX 

Uyy + 2(^11 — Kl2)u^y\ 

+ \[n{dref + piid^ef] 

+ a±daUaidrO 

+ \[Pn{.drUzf + Kii/A{d^Uzf] (33) 
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where a, /3 = a;, y. By comparing Eqn l33l with Eqn l271 
we can see that Kn A + 2/i, K12 — ?> A, so all the 
discussions in the isotropic case can be used here after 
the replacements. Fig. 2 follows after these replacements. 
Namely, only the longitudinal component in the xy plane 
is coupled to the 6 mode, while the transverse mode in 
the xy plane with velocity v^^y = {Kn — Ki2)/2pn is 
decoupled. Obviously the transverse mode along z di- 
rection Uz mode with the velocity vl^ = K44/4pn is also 
decoupled. Note that the two transverse modes have dif- 
ferent velocities. It was found that vip ~ 455m/ s,Vtz ~ 
255m/ s, vtxy ~ 225m/ s when q is along the xy planed. 

Along any general direction g, strictly speaking, one 
can not even define longitudinal and transverse modes, so 
the general action Eqnl^ should be usedSS. Despite the 
much involved 4x4 matrix diagonization in Ux,Uy,Uz, 0, 
we expect the qualitative physics is still described by 
Fig.2. 

In principle, inelastic neutron scattering experiments'^^ 
or acoustic attenuation experiments^ can be used to de- 
tect the predicted the low energy excitation spectra in 
the SS shown in Fig.2. 



VIII. TOPOLOGICAL ELEMENTARY 
EXCITATIONS IN SS: VORTEX LOOPS AND 
VORTEX 

In the last section, we studied the low energy excita- 
tions shown in the Fig.l by neglecting the topological 
vortex loop. Here, we will study how the vortex loop in- 
teraction in SS differ from that in the SF. For simplicity, 
in the following, we only focus on the isotropic case. The 
formulations can be generalized to the hep case straight- 
forwardly. We can perform a duality transformation on 
Eqnl24]to the vortex loop representation: 

JOv = 7rj^{^ti.v\crdua\a - adaUaS^rT + ^STTa^i^j^^ (34) 

where /i, u, A, a stand for space and time, but a, /3 stand 
for the space components only, Kq = k, Ka = Ps and 
a^jy = — fljyp is an anti-symmetric tensor gauge field and 
jj^jj — -^(-fivXadxduO is the anti-symmetric tensor vortex 
loop current due to the topological phase winding in 6. 

Eqn[M] has the gauge invariance a^i, — >■ a^v -\- d^Xf ~ 
di^Xfj. where Xn is any 4-component field'*^. It is the most 
convenient to choose the Coulomb gauge datta^ = to 
get rid of the longitudinal component, then the transverse 
component is at = iSap-yQa^p-y/l- It can be shown that 
|atp = 2|aa^p. Then Eqnl31is simplified to: 

£y = i[p„a;,^j + (A + 2^ + a^/K)g^]|u;(g,w„)P 

+ liqV^ + ^l/ps)\atf + —q^\aoc.f 
2 Ps 

- aq^/KUi{-q,-uJn)at{q,i^n) 

+ i27rjo^aoa + i2Trj'^paap (35) 



where the transverse phonon mode ut was omitted, be- 
cause it stays the same as in the NS as shown in Eqnl?7l 
It is easy to see that only at has the dynamics, while 
aoa is static. This is expected, because although a^j/ has 
6 non- vanishing components, only the transverse com- 
ponent at has the dynamics which leads to the original 
gapless superfluid mode = v'^q'^. Eqn[551 shows that 
the coupling is between the longitudinal phonon mode 
ui and the transverse gauge mode at- The vortex loop 
density is Jq^ = -^eap-ydpd^O and the vortex current is 
fap = ^ea/37[9o,9T,]6'. Integrating out the aoa, we get 
the vortex loop density-density interaction: 

TiPs I dr [ dxdyj^^{x,T)--^—j^^{y,T) (36) 

Jo J F ~ y\ 

Namely, the vortex loop density- density interaction in 
SS stays as l/r which is the same as that in NS ! There- 
fore, a single vortex loop energy and the critical transi- 
tion temperature T^dxy in Fig.l is solely determined by 
the superfluid density ps independent of any other pa- 
rameters in Eqn[27l except that the vortex core of the 
vortex loop is much larger than that in a superfluidi^. 

In a cylindrical geometry used in the torsional oscil- 
lator experiment, the vortex loops will become straight 
vortex line alone the rotational axis. In the SF phase, 

a single vortex costs a lot of energy E^^ — \ In-^ 
where m is the mass of He atom, R is the system size 
and ^sF a is the core size of the vortex. This energy 
determines the critical velocity in SF vf^ > 30cm/s. 
Because the long distance behavior of SS is more or less 
the same as SF, we can estimate its single vortex energy 

E^^ — where pf^ is the global superfluid 

density inside the SS. We expect the core size of a super- 
solid vortex ^l/K^l/kr^a^S_sF- So inside 
the SS vortex core, we should also see the lattice struc- 
ture of r^. This is similar to the phenomenon that DW 
ordered states were detected in the vortex core of high 
temperature superconductors^Si^. In fact, because i^{x) 
stands for vacancies or interstitials, we expect that S^ss 
should be of the order of the average spacing between the 
interstitials or vacancies in the SS. It is interesting to see 
if neutron or light scattering experiments can test this 
prediction. Compared to E^^ , there are two reductions, 
one is the superfluid density, another is the increase of 
the vortex core size ^55 » ^sf- These two factors con- 
tribute to the very low critical velocity vf^ ~ 30pm/s. 
Of course, the reduction from the increase of the vortex 
core is negligible because of the logarithmic dependence. 

Integrating out the aap , we get the vortex loop current- 
current interaction: 

'^'^'^fapi-q^ ~^n)Dap,js{q,UJn)jl!^s{q^'^n) (37) 
where Dap,'ys{q,UJn) = {Sa-yS/SS - Sp-^SaS - ^^^a^ - 

^Sf}s + ^Sa5 + ^S^s)Dt{q,iUn) whcrc Dtiq,ujn) is 
the at propagator. Defining ADt{q,uJn) = Df^ {q,uJn) — 
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Df^{q,ujn) as the difference between the at propagator For simphcity, we just give the expression for the equal 
in the SS and the SF, then from Eqn l35[ we can get: time 

? ? 4 

^jj^ ^ (38) 



ADt{x 



0) 



1 



(39) 



Namely, the vortex current-current interaction in SS is 
stronger than that in the SF with the same parameters 



IX. X-RAY SCATTERING FROM THE SS 

Let's look at the prediction of our theory on X-ray 
scattering from the SS. For a lattice with j = 1, • ■ ■ i 
basis located at dj , the geometrical structure factor at the 
reciprocal lattice vector K is S{K) — J2^=i fj{K)e^^''^^ 
where fj is the atomic structure factor of the basis at dj . 
The X-ray scattering amplitude I{K) ^ \S{K)\'^ . Again, 
we discuss SS-v and SS-i respectively. 

A. X-ray scattering from the SS-v 

In this case, because the superfluid density wave simply 
sits on the n lattice, so the X-ray scattering is very similar 
to that from NS at mean field level. However, as shown 
in C, quantum and thermal fluctuations will still make 
the X-ray scattering from the SS-v different from that 
from the NS. 



B. X-ray scattering from the SS-i 

In this shown in section VI, the superfluid den- 

sity wave is shifted from n lattice, so the X-ray scattering 
is different from that from NS even at mean field level. 
For simplicity, we first take the sc lattice to explain the 
main points, then list the X-ray scattering from all the 
other lattices classified in section VI. 

(a) Simple Cubic lattice. For the SS in the sc lattice, 
as shown in (c2) of the last section, the local superfluid 
density attains its maximum at the dual lattice points 
of the sc lattice. Then di = 0,Zi = §(?+ J+ k),K = 
^{nii+n2j+n3k), then taking the ratio of the geometric 

structure factor of SS over that of the NS Sss [K) /Sns — 
1 + /(-l)"i+"2+"3 where / - - a^. It is 1 + / for 
even K and 1 — / for odd K. 



(b) Triangular lattice, di = 0, Zi = \{ai + 02), ^2 = 

|(ai + 0,2), K = ni&i + 712^2, then taking the ratio of 
the geometric structure factor of SS over that of the NS 
Sss{K)/Sns = l + 2/cos^(ni-hn2). This result could 
be relevant to possible 2d excitonic superfluid in electron- 
hole bilayer system to be briefly mentioned in section X. 

(c) hep lattice, hep lattice is not a Bravais lattice. In 
NS, di =0,^2 = i(ai+a2)+a3/2. A' = ni6i+rt2?2+"-3&3, 
then Sns{K) = 1 + e'2'^(^^4^+T^). In the SS, there are 
2 more additional basis at li — |(ai -I- 0,2) + 123/4, /2 = 
803/4, then Sss{K) = Sns{K) + /e'2.(^^^4^+^) + 
/e'¥"3 ^ Sns{K) + fe-'i^^S*^s(K). 

(d) bcc lattice. We think bee lattice as a sc lattice plus 
a basis, di = 0^d2 = ^{i + j + k),K = ^{nii + n2j + 
nsk), Then SNsiK) = 1 + (-l)»i+«2+n3 ^Y^i^h is 2 for 
even K and for odd K. In the SS, there are 3 more 
additional basis at = |(i + j),/2 = |(i + fc),Z3 = + 

fc), then5s5(^) = 5Ar5(i?) + /[(-l)"^+"^ + (-l)"^+"-' + 

(e) fee lattice. We think fee lattice as a sc lattice plus 
4 basis located at di ^ 0,d2 = ^{i+j), d^ = ^4 — 
^ij + k),K = ^{mi + n2j + nak), then Sns{K) = 
l-H[(-l)"i+""+(-l)"i+"-^ + (-l)""+"^']. In the SS, there 
is one more additional basis located at = |(i + j + fc), 
then SssiK) = SNsiK) + /(-l)"i+"2+n3. it is easy to 
see that in bcc and fee lattices, we need simply exchange 
d vectors for the NS and the / vectors for SFDW. 

Wc conclude that the elastic X-ray scattering intensity 
from the SS-i has an additional modulation over that 
of the NS. The modulation amplitude is proportional to 
the maxima of the superfluid density ~ which 

is the same as the NCRI observed in the PSU's torsional 
oscillator experiments. 



C. Debye- Waller factor in the X-ray scattering 
from the SS-v and SS-i 

It is known that due to zero-point quantum motion in 
any NS at very low temperature, the X-ray scattering 
amplitude I{G) will be diminished by a Debye- Waller 
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(DW) factor ~ e~3<^^^"°^ where Ua is the lattice phonon 
modes in EqnlJH In Eqni24l if the couphng between 
the u and were absent, then the DW factor in the SS 
would be the same as that in the NS. By taking the ratio 
Iss{G)/ Ins{G) at a given reciprocal lattice vector G, 
then this DW factor drops out. However, due to this 
coupling, the (u^) in SS is different than that in NS, so 
the DW factor will not drop out in the ratio. In this 
subsection, we will calculate this ratio and see how to 
take care of this factor when comparing with the X-ray 
scattering data. 

As identified below EqulHl the density order param- 
eter at the reciprocal lattice vector G is pg{x,T) = 
e^G-u(x,T)^ then (pg(f,T)) = e-5'^-G,(«,u,)_ r^j^g Debye- 
Waller factor: 



where the phonon-phonon correlation function is: 

{uiUj) = {uiui)qiqj + {utUt){5i-j - qigj) (41) 

where qicjj = . 

Then substituting Eqn|4T]into Eqn|40] leads to: 

a{G) = IssiG)/lNsiG) = e-3G=[K(^.^)>ss-(uf (5,r))«s] 

(42) 

where the transverse mode drops out, because it stays 
the same in the SS and in the NS. 

Defining {Au'^)i{q,iLOn) = (|u;(g, ia;„)p)55 - 

{\ui{q,iuJn)\^)NS, {Au'^)l{cj) = Y.i^S^U^W^i^n) 



and {Au 



uj{x,T)) 



I(G) = \{pMx,T) 



^ — GiGj {ui{x,T)uj {x.t)) 



(40) is easy to see: 



SS - {uj{x,t))ns = 
= J^iAu'Uq), it 



{Au^)i 



d?q 1 



E 



-a q^n 



(27r)3 p ^ {[^l + p.g2)(p„^2 + (A + 2^)q2) + a\^ul\{p,,u:l + (A + 2p)q^\ 



(43) 



r 



Obviously, (Am^); < 0, namely, the longitudinal vibra- 
tion amplitude in SS is smaller that that in NS. The 
a{G){T 0) = g-|G^(A«^), ^ I rpj^ig expected, be- 
cause the SS state is the ground state at T < Tss, so the 
longitudinal vibration amplitude should be reduced com- 
pared to the corresponding NS with the same parameters 

Pn,\lJ- 

After evaluating the frequency summation in Ean l43[ 
we get: 



(Pq 1 coth/3i)+q/2 coth^Upg/2 



(27r)3 p„ 2v+q 2vpq 
v'^-v'^ coth(3v+q/2 coth;3w_(j/2 



-)(- 



2v+q 



2v^q 



At T = 0, the above equation simplifies to: 

d^q I , 1 1 



{Au%iT = 0) 



(27r)3 pn 2v+q 2vpq 



-){-. 



v'^ — v'^'''2v+q 2v^q' 

_ {V+ +V- ~Vp- Vs) A^ 

1 A^ 

npn (W+ +Vp+ Vs)iv+ + V-)Vp STT^Pn 



< 0(45) 



where A ^ 1/a is the ultra-violet cutoff and we have used 
the fact v+ + > Vp + Vs- 



By subtracting EqnHS] from EgnEl we get 

{Au^)i{T) - iAu^)i{T = 0) = 
{v+ - Vp){v+ + Vp) - {vs - w-)(w+ + V-) {ksTY 



(u+ -|- V-)v+V-Vp 



12 Pn 



> (H6) 



Namely, the difference in the ratio will decrease as as 
the temperature increases. Of course, when T approaches 
Tss from below, the difference vanishes, the a{G) will 
approach 1 from above, the SS turns into a NS. 



D. Density-density correlations 

The density-density correlation function in the SS is: 

{p^ix,t)p*^{x',t')) ^ e-iG,G,{iu,ix.t)-u,ix'.t')){u,{x.t)^u,{S',t' 

(47) 

where t is the real time. 

For simplicity, we only evaluate the equal-time corre- 
lator {pg{x,t)p*^{x',t)) = {p^{x,T)p*g{x',T)) where r is 
the imaginary time. It is instructive to compare the den- 
sity order in SS with that in a NS by looking at the ratio 
of the density correlation function in the SS over the NS: 

apix^x') = {PaPQ)ss/{PGP*G)NS = 

It is easy to find that 
d^q 



(48) 



ADp{x-^) 



(27r) 



-(2-e 



iq-{x—x') ^—iq-{x—x')\( \ ^ ,2 



(49) 
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where {/^v?)i{q) is defined above Eqnl43l and is the inte- 
grand in EqnHH 

At r = 0, the above equation can be simphfied to 



1 



1 



(50) 

So we conclude that ap{x—x') < 1, namely, the density 
order in SS is weaker than the NS with the corresponding 
parameters p„, A, /i. This is expected because the density 
order in the SS is weakened by the presence of moving 
vacancies. 

Unfortunately, so far, the X-ray scattering data is lim- 
ited to high temperature T > 0.8K > Tsa^. X-ray 
scattering experiments on lower temperature T < Tss 
are being performed to test these predictions. 



X. 



SPECIFIC HEAT IN THE SS 



It is well known that at low T, the specific heat in 
the NS is = + CP/ + C^an where = 

"^j^kBiJj^Y is from the longitudinal phonon mode and 

Gfp^ = 2 X ^fcsdt^)^ is from the two transverse 
phonon modes, while C„a„ is from the vacancy contri- 
bution. Cvan was calculated in-^ by assuming 3 different 
kinds of models for the vacancies. So far, there is no 
consistency between the calculated C^an and the experi- 
mentally measured one^'^''^^. The specific heat in the SF 
C^P = ^kei^f is due to the SF mode 0. In this 
subsection, we focus on the specific heat inside the SS. 
From Eanl27l we can find the specific heat in the SS: 



SS 



2^ 
15 



ksi'^r+'-^ksi'f^f + C^^ (51) 

nv4- 15 nv- 



where C*p stands for the contributions from teh trans- 
verse phonons which are the same as those in the NS ( 
Fig.7 ). 

It was argued inlS,, the critical regime of finite temper- 
ature NS to SS transition in Fig.l is much narrower than 
the that of SF to the NL transition, so there should be a 
jump in the specific heat at T = T^S* ( Fig.7 ). EanlSTl 
shows that at T < Tss, the specific heat still takes ~ 
behavior and is dominated by the w_ mode in Fig. 6. 

From Eqn l511 it is easy to evaluate the entropy inside 
the SS: 



SssiTss)-^-^kB{ ^ 



''^^''n^ + ^) + S^^ (52) 



where S'*^ is the entropy due to the 2 transverse modes. 
The excess entropy due to the vacancy condensation 

is: 



45 



ksTss ^ 1 



1 



(53) 
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FIG. 7: The specific heat of a supersolid 

Obviously, A5 > is due to the lower branch < 
vip. At a = 0, the above equation reduces to AS = 
^^ksi '^^j^^^ )^ which is simply the vacancy condensa- 
tion into SF. Using the molar volume vq ~ 20cm^ /mole 
of solid '*-ffe and Tss ^ lOOmK, we can estimate the AS" 
per mole is ~ 10~^R where R is th gas constant. This es- 
timate is 3 orders magnitude smaller than that in^ where 
the SS state was taken simply as the boson condensation 
of non-interacting vacancies. Our estimate is indeed con- 
sistent with recent experiment on specific heat^. The 
linear ^ T behaviour found in heat capacity experiment 
can only be due to disorder or ^He impurities. 



XI. ELEMENTARY EXCITATIONS IN FULDE - 
FERRELL-LARKIN-OVCHINNIKOV (FFLO) 
STATE OF SUPERCONDUCTORS 

The QGL theories constructed for the supersolid 
in the previous sections may also be used to study 
the elementary excitations in the inhomogeneous 
Fulde - Ferrell-Larkin-Ovchinnikov (FFLO) state of 
superconductors^"— . As argued in the following, the 
FFLO state can be taken as a fermionic supersolid. When 
the number of spin up electron is equal to the number of 
down spin electron n-^ = n^, the pairing is at fc = only. 
If there is a mismatch Sn = — n^, then pairing will 
shift to a non-zero momentum go — kp-f ~ kp^. By us- 
ing the GL theory near the normal to the FFLO state 
( Fig. 8 ), at a mean field level, the authors in^ found 
the most favorable lattice structures of the FFLO state 
is the stripe state (LO state) in large number of param- 
eter regimes. The FFLO state maybe considered as a 
weak coupling ( or fermionic ) analog of the ( bosonic ) 
supersolid. Indeed, the Fig. 8 is similar to the Fig. 2 after 
identifying the chemical potential difference with the 
pressure p. It is important to point out a few important 
differences between the GL for the bosons developed in 
this paper and the GL for the FFLO state used in^: (1) 
In the former, due to the lack of the Z2 symmetry, there 
is a cubic term in Eqnjl] But in the latter, due to the 
presence of the Z2 symmetry, there is no such a cubic 
term. (2) In the latter, the transition from the normal 
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FIG. 8: The phase diagram of FFLO state. Sfj, is the difference 
of the two chemical potentials of two species of fermions, T is 
the temperature, Aq is the energy gap at the balanced case 
6n = 0. 



state to the FFLO state is a classical Lifshitz type tran- 
sition, so there is no zero momentum pairing as shown 
in^^ and in Eqn l54[ in contrast to the second equation 
in EgnllQI where there is always a zero momentum BEC 
condensation. (3) In the latter, there are large number 
of unpaired fermions which are the normal component of 
the system. Their distributions also take the same lattice 
structure as the FFLO state. They play a similar role as 
the normal lattice component in the first equation in 1191 
but they form extended Bloch waves in the underlying 
FFLO lattice and are integrated out in the GL theory 
irtSL. 

So far all the analysis in a FFLO stated"— are only 
at a mean field level. In this appendix, we study the ele- 
mentary excitations above the mean field solutions. Just 
from symmetry breaking point of views, the FFLO state 
breaks both the translational symmetry and the U{1) 
symmetry, so it also has both the diagonal order and off- 
diagonal order, so the low energy modes should contain 
both the lattice phonon mode and the superfluid phonon 
mode. The pairing order parameter can be written as: 



FFLO 



{x, r) = Vo X! ' 



(54) 



where G — qq, the tpo{x,T) = Ae'^'-^''^' is the pairing 
order parameter, u{x, r) are the lattice phonon modes. 
Note the absence of the zero momentum pairing. For 
a charged system such as a electron system, due to the 
Higgs mechanism, the Goldstone mode d{x, r) will be just 
eaten by the gauge field. However, for a neutral system 
such as the pairing between two species of fermions with 
unequal populations in cold atom systems across a Fesh- 
bach resonance^^, the Goldstone mode 6{x, r) stays. The 
coupling between the two phonon modes in Eqnl54] are 
also described by a equation similar to EgnlTTl The only 
difference is that the lattice structure is a LO state in- 
stead of a isotropic solid. After taking this difference into 
account, the elementary excitations inside the FF state 
are similar to those in the Fig. 6a and the corresponding 
spectral weights can be worked out similarly. The exper- 
imental signatures of the elementary excitations can also 



be worked out similarly. 

There exist other kinds of fermionic supersolids in ad- 
dition to the FFLO state. In some fermionic systems, 
it is easy to see the coexistence of CDW and Super- 
conductivity ( SO), so " fermionic supersolid " phases 
are common. For instance, a quasi-two-dimensional sys- 
tem NbSe2 has a transition to an incommensurate CDW 
phase at some high temperature Tcdw and then a tran- 
sition to a phase with coexisting CDW and SC order 
at a lower temperature Tsc- The CDW is a pairing 
in particle-hole channel at 2fci?, its order parameter is 
"ipCDW — {Ca{k)ca{k + Q)) where Q is the ordering wave 
vector of the CDW. The SC is a pairing in particle- 
particle channel also across the Fermi surface. Its order 
parameter is ipse = (c|(^)cj(— fc)). Both order parame- 
ters are composite order parameters. Different parts of 
Fermi surface can do the two jobs separately ( see, for 
example)^ ). In contrast to the bosonic SS where there 
is a density operator n and a complex order parameter "0, 
both order parameters ipc dw and ^psc are complex order 
parameters. We can see that although the formation of 
a supersolid in a bosonic system has completely different 
mechanism than in fermionic case, from symmetric point 
of view, the QGL theories constructed developed in pre- 
vious sections could also be used to describe the interplay 
between the two complex order parameters, the proper- 
ties of the "fermionic supersolid" (FSS) and the transi- 
tion from the FSS to the CDW. However, the fermionic 
excitations near the nodes maybe also be important and 
need to be taken into account"*^. 



XII. CONCLUSIONS 

In this paper, assuming there are quantum fluctuations 
generated vacancies or interstitial at zero temperature^, 
we constructed a two component QGL theory to map out 
a possible global phase diagram, analyze carefully the 
conditions for the existence of the supersolid and study 
all the phases and phase transitions in a unified view. 
The only new parameters introduced in the GL theory 
in this paper is the coupling g and v between the n sec- 
tor ( or normal solid part ) and the ip sector ( or the 
superfluid part ) in EqnlSl We investigated the SS state 
from both the SF and the NS side and found completely 
consistent description of the properties of the SS state. 
Starting from the SF side with increasing the pressure, 
we developed the theory basing on the two facts (1) there 
is a roton minimum in the superfluid state (2) the insta- 
bility to solid formation is driven by the gap diminishing 
at the roton minimum. By increasing the pressure from 
the superfluid side, there are two possible scenarios (1) 
the SF to the C-NS transition in Fig.l was described by 
the QGL action EqnlT^lfirst derived in this paper (2) the 
SF to the SS transition in Fig. 2 is a simultaneous combi- 
nation of the SFDW transition in the ip sector driven by 
the roton condensation at ko = kr and the NS transition 
in the n sector driven by the divergence in the structure 
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function ko = kn = kr- The superfluid becomes a SS 
at lower temperature and a NS at higher temperature ( 
Fig. 2 ). Then we also approached the SS state from the 
NS side. Depending on the sign and strength of the cou- 
pling g between the solid and superfluid, we also found 
two possible scenarios: (1) If \g\ is sufficiently small ( 
Fig. 4 ), then the resulting solid at T = 0,pci < p < Pc2 
is a commensurate normal solid ( C-NS ). The SS state 
does not exist as a ground state. However, it may still 
exist as a metastable state. The QGL action to describe 
this SF to NS transition in Fig.l was developed in Sec. II. 
(2) If \g\ is sufficiently large ( Fig. 4 ), the resulting solid 
at T = 0,Pci < p < Pc2 is an incommensurate solid with 
zero point quantum fluctuations generated vacancies if it 
is negative and interstitials if it is positive ( Fig. 4 ). The 
condensation of the vacancies or interstitials lead to the 
formation of the SS-v and SS-i respectively. The SS state 
has lower energy than the NS state at T = 0. The Tss-v 
( Tss-i ) is an effective measure of the strength of the 
interaction g in the SS-v ( SS-i ) supersolid. There is no 
particle-hole symmetry relating Tss-v to Tss-i- Many 
physical consequences came out of this single parameter 
g. Our results on supersolid should be independent of 
many microscopic details and universal. 

Just like the SF is a uniform two-component phase con- 
sisting of superfluid and normal component at any finite 
temperature, the SS state is a uniform two-component 
phase consisting of a superfluid density wave (SFDW) 
and a normal solid component even at zero temperature. 
The SFDW in the SS-v coincides with the underlying 
normal solid. While the SFDW in the SS-i state is just a 
dual lattice to the underlying normal solid. This impor- 
tant fact leads to one of the key predictions in this paper: 
the X-ray scattering intensity from the SS-v is similar to 
that of NS at mean field level, while the X-ray scattering 
intensity from the SS-i ought to have an additional mod- 
ulation over that of the NS. The modulation amplitude 
is proportional to the Non-Classical Rotational-Inertial 
(NCRI) observed in the torsional oscillator experiments. 
However, the X-ray scatterings from SS-v and SS-i will 
be modified by the Debye- Waller factor calculated in sec. 
IV. The NS-v ( NS-i ) to SS-v (SS-i ) transition is de- 
scribed by a 3d XY model with much narrower critical 
regime (Fig. 4). 

We also studied the zero temperature quantum phase 
transition from the SS to the NS driven by the pressure 
near the upper critical pressure p = Pc2 in Fig. 5. We 
found that the coupling to the quantum fluctuation of 
the underlying lattice is irrelevant, so the transition stays 
the same universality class as the superfluid to Mott in- 
sulator transition in a 3 dimensional rigid optical lattice. 
The flnite temperature transition from the SS to the NS 
in Fig.l was studied previously in^^ and in'^'^ in differ- 
ent contexts. It was found that the coupling to classical 
elastic degree of freedoms will not change the universal- 
ity class of the 3D XY transition. However, we found 
that the coupling to quantum lattice phonons is very im- 
portant inside the SS and leads to two "supersolidon" 



modes uj± = v±q { one upper branch and one lower 
branch ) shown in Fig. 2. Their corresponding spectral 
weights are also worked out. The transverse modes in 
the SS stays the same as those in the NS. Detecting the 
two supersolidon modes with the corresponding spectral 
weights, especially, the lower branch aj_ mode by neutron 
scatteringiZ. or acoustic wave attention experiments'^^ is 
a smoking gun experiment to prove or disprove the ex- 
istence of the SS in helium '^He. The tu- is estimated 
to considerably lower than the sound speed in the super- 
fluid. Then we calculated the experimental signature of 
the two supersolidon modes. We found that the longitu- 
dinal vibration in the SS is smaller than that in the NS 
( with the same corresponding solid parameters ) , so the 
Debye- Waller factor at a given reciprocal lattice vector 
is larger than that in the NS. The density-density corre- 
lation function in the SS is weaker than that in the NS. 
By going the to the dual vortex loop representation, we 
found the vortex loop density-density interaction in SS 
stays the same as that in the SF ( with the same cor- 
responding superfluid parameters ), so the vortex loop 
energy and the corresponding SS to NS transition tem- 
perature is solely determined by the superfluid density 
and independent of any other parameters. The vortex 
current-current interaction is stronger than that in the 
SF. The specific heat in the SS is still given by the sum 
from the transverse phonons and the two supersolidon 
modes and still shows ^ behaviors. The supersoli- 
don part is dominated by the lower branch. The NCRI 
is only weakly anisotropic in the SS phase for hep lat- 
tice. In principle, all these predictions can be tested by 
experimental techniques such as X-ray scattering, neu- 
tron scattering, acoustic wave attenuation s '^^''^^ and heat 
capacity. 

It may be instructive to make some analogy of Fig.l 
and Fig. 2 at T < Tsf to Type-I and type-II supercon- 
ductors with the pressure p playing the role of the mag- 
netic field H: Fig.l is similar to Type-I superconductor 
with SF identified as the Messiner state, the NS as the 
normal state, the critical pressure pc identified as the 
critical magnetic field He- Fig. 2 is similar to Type II 
superconductor with SF identified as the Messiner state, 
the SS as the mixed vortex lattice state which also breaks 
both translational order and the global U{1) symmetry, 
the NS as the normal state, the lower and upper critical 
pressures pd and Pc2 identified as the lower and upper 
critical magnetic fields Hd and Hc2- In superconductors, 
it is the K = A/^ to determine Type I and Type II and 
if the vortex lattice is a stable intermediate state or not 
as the magnetic field is increased. In Helium 4, it is the 
sign and strength of the coupling constant g in Eqnl3] 
to determine the Fig.l and Fig. 2 and if the SS is a sta- 
ble intermediate state or not as the pressure is increased. 
So the pressure p and the coupling g in the formation of 
SS-v play the role of the magnetic field H and k in the 
formation of the mixed state of superconductors. Note 
that in superconductors, H and n are two independent 
parameters, in '^He, p and g are also two independent 
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parameters. 

The GL theory developed in this paper put the com- 
peting orders of superfluid and sohd in the unified frame- 
work. We suggest that even supersohd may not be re- 
ahzed in '^He system, it has its own intrinsic, scien- 
tific interests and may be reahzed in other continuous 
bosonic and fermionic systems. For example, in symmet- 
ric electron- hole bilayer systemsiSi"—, it was shown in*^^ 
that it is quite possible that there may a narrow window 
of ESS where both order parameters are non-vanishing 
(tp) ^ 0, (ng) 7^ intervening between the ESF and 
two weakly coupled Winger crystal. A similar GL theory 
can be used to study the elementary excitations inside a 
FFLO state and some fermionic systems with coexistence 
of CDW and superconductivity. 
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Appendix A: Discussions on a tight-binding toy 
Supersolid ground state wavefunction 

The Ginzburg-Landau theory constructed in the main 
text is based on order parameters and symmetries. It 
should hold irrespective many microscopic details such 
as what is the mechanism responsible for the formation 
of the supersolid. Despite there are many microscopic 
calculations for ^i/e, constructing a microscopic theory 
for supersolid is very difficult. In this appendix, I will 
discuss a well known toy SS wavefunction and clarify a 
few concepts related to global phase-number uncertainty 
relation and the role of vacancies or interstitials in the for- 
mation of SS. We also clarify the physical meaning of the 
order parameters n in EqnlT]and ip in Eqn. [21 However, 
because the toy wavefunction may miss some important 
physics in bulk ^He systems. For example, due to the 
very peculiar potential well in the solid '^He which has a 
local shallow maximum at the lattice site, the tight bind- 
ing model maybe crude, so the discussion is intuitive and 
instructive. Eqn lAll is also a ground state wavefunction, 
so does not include the lattice phonon interaction. 

The toy wavefunction of a supersolid takes the BCS 
like form 

N 

\SS)=l[{u + vbl)\0} (Al) 

i=l 



where N is the number of sites, m ^ and 
If setting u — 0, the state reduces to a commensurate 
solid without any vacancies \CS) = Hi^i^IlO)- The 
commensurate solid (CS) is an exact eigenstate of the bo- 
son number operator N}, = X]i3=i '^j with the eigenvalue 
Nb = N, so has no chance to become phase ordered. 
Adding a superfluid component to the CS leads to the 
SS in Eqn IAll If setting v = 0, then the state reduces to 
the vacuum state |0). 

If setting u — cos|,w = sin|e"^, then Eqn IAll be- 
comes: 

'^0 
\SS,cl^}=l[{cos-+sin-e^Hlm (A2) 

where 9 tt. 

By construction, the state has the translational order 
with the average boson density (ui) = \v\'^ = sin^0/2, 
so the average vacancy density is jitp = cos^6'/2. The 
total number of bosons Nt — Nsin^ 9/2, the number of 
vacancies = N cos^ 9 /2. It is easy to see that \SS) 
also has the ODLRO with (hi) — u*v — ^siii9e'^^, so 
\SS) is indeed a supersolid state. The angle 9 controls 
the magnitude, while the phase controls the phase of 
the condensation. Defining bk=o = J2iLi which 
satisfy the boson commutation relation [fco7&o] — 1, the 
boson operator at zero momentum is no = b\=o^k=o, 
the total number of bosons at the zero momentum 
state is Nq ^ {SS\nfi\SS) = f sin^e* = NbCOs'^9/2 < 
Nb = Nsm^9/2 < N. Similarly iVo = (S'S'|no|S'5) = 
f sin^ei = Ny 3111^9/2 < ^ Ncos'^9/2 < N. So 
the condensate is less than both the boson density and 
the vacancy density. At integer filling n = 1, the 
non-interacting BEC state \SF) = -^(6|,^q)^|0), then 

iVo = {SF\no\SF) ^ Nb ^ N. Obviously, this non- 
interacting BEC state is not included in the family in 
the Eqn El 

A supersolid state \SS, Nb) with Nb bosons is given by: 
\SS,Nb)= J^'^ ^e-^''^*\SS,cl,) (A3) 

where the total boson number Nb and the global phase 
(j) are two Hermitian conjugate variables satisfying the 
commutation relation: [Nb, (j)] — ih. It leads to the un- 
certainty relation ANbAip > 1. 

It is easy to see (Nb) = N\v\'^ = Nsm^9/2,ANb = 
V(A^^) - {my = VN\uv\ = VN^\sm9\,ANb/{Nb) = 
--l=\u/v\ — -^|cot0/2|. E 9 ^ TT, the absolute boson 

number fluctuation ANb ^ Vn is quite large, so A9 
could be quite small, so one can get a phase ordered 
state. On the other hand, the relative boson number 
fluctuation ANb/{Nb) ^ is quite small, so one can 
still measure the average boson number accurately. The 
first quantization form of the Eqn lA3l can be derived by 
the same method used 
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Because in the SS state, there is a global phase order- 
ing in (j), so its conjugate variable is the total number of 
particles as shown in Eqn lA3l The local tunneling or 
exchanging processes stressed in** may not cause the total 
number fluctuations, therefore may not cause the global 
phase ordering leading to the supersolid phase. Eqn IAll 
implies that the vacancies in an incommensurate solid 
could lead to the formation of a supersolid. If this indeed 
happens, in the GL theory constructed for SS-v in the 
main text, the lattice sites are represented by n(a;), while 
the vacancies are represented hy tp. Of course, the trial 
ground wavefunction lAll does not include fluctuations, 
so no information on the lattice phonons and superfluid 
phonons, therefore no supersolidons in the Fig. 6a. A toy 
wavefunction for SS-i was not written down so far, be- 
cause the interstitials are moving between lattice sites, 
so there is no tight binding limit. The interstitial case is 
described by the extended Bloch state in Section IV-B. 
In the SS-i, the bosons are represented by n{x), while the 
interstitials are represented by "0. 



Appendix B: Comparisons with Supersolids on 
lattices 

The extended boson Hubbard model ( EBHM ) with 
various kinds of interactions, on all kinds of lattices and 
at different filling factors is described by the following 
Hamiltoniar^: 

H = -t ^(6^6, + /i.e.) + ^ 51 

where Ui = b^bi is the boson density, t is the near- 
est neighbor hopping amplitude. U,Vi,V2 are onsitc, 
nearest neighbor (nn) and next nearest neighbor (nnn) 
interactions respectively, the • • • may include further 
neighbor interactions and possible ring-exchange inter- 
actions. A supersolid in Ean lBll is defined as the simul- 
taneous orderings of ferromagnet in the XY component 
( namely, (bi) ^ ) and CDW in the Z component. 
In^, we studied all the possible phases and phase tran- 
sitions in the EBHM in bipartite lattices such as a hon- 
eycomb and square lattice near half filling. We show 
that there are two consecutive transitions at zero tem- 
perature driven by the chemical potential: in the Ising 
limit, a Commensurate-Charge Density Wave (CDW) at 
half filling to a narrow window of CDW supersolid, then 
to an Incommensurate-CDW ; in the easy-plane limit, 
a Commensurate- Valence Bond Solid (VBS) at half fill- 
ing to a narrow window of VBS supersolid, then to an 
Incommensurate-VBS. The first transition is second or- 
der in the same universality class as the Mott to insula- 
tor transition , therefore has the exact critical exponents 
z = 2,1/ = 1/2, 7] = with logarithmic corrections, while 



the second one is first order. Liu and Fisher^* also stud- 
ied the C-CDW to the CDW-SS transition and concluded 
that z = 1 in contrast to z = 2 achieved in^. We found 
that the phase diagram in the Ising limit is similar to the 
reentrant " superfluid" in a narrow region of coverages in 
the second layer of '^He adsorbed on graphite detected 
by Crowell and Reppy's torsional oscillator experiment 
in 1993^"^. Indeed, the data in the torsional oscillator 
experiment in^ do not show the characteristic form for 
a 2d *ife superfluid film, instead it resembles that 
characteristic of a supersolid in terms of the gradual on- 
set temperature of the NCRI, the unusual temperature 
dependence of Tss on the coverage. Of course, both ex- 
periments may due to phase separations instead of the SS 
phase, very recently, the author studied various kinds od 
supersolids in frustrated lattices such as triangular and 
kagome lattices^. 

By comparing Eqnsll] I2|3I in the continuum with the 
Egn lBll on the lattice, one can see the crucial difference 
between a lattice system with a continuous system: Eqn 
[1] stands for a spontaneously formed lattice with phonon 
excitations. While Eqn IBll is a fixed external lattice with 
no such phonon excitations. This crucial difference is 
responsible for the differenceSS in the elementary exci- 
tation spectra in the two kinds of supersolid shown in 
Fig. 6. So the reentrant lattice SS discussed i n^^iSOiSi jg 
different from the bulk '^He SS state discussed in this 
paper, although both kinds of supersolids share some in- 
teresting common properties. In both cases, the SF to 
SS transition is driven by the closing of the gap of the ro- 
ton minimum, so the transition could be first or second 
order at mean field level. However, due to the fluctu- 
ation of the underlying superfluid order, the transition 
turns out to be the flrst order as conflrmed by the QMC 
simulations. On the lattice, the manifold of the roton 
minimum consists of discrete points due to the lattice 
symmetry. However, in '^He, as shown in sections HI 
and IV, the manifold of the roton minimum is a continu- 
ous surface, so the transition must be first order, order. 
In the former, there is a periodic substrate or spacer po- 
tential which breaks translational symmetries at the very 
beginning. The filling factor is controlled by an external 
chemical potential. At integer filling factors or half filling 
factor, there are particle-hole symmetry (PH) for excita- 
tions which ensures the number of particles is equal to 
that of vacancies and adding interstitials is equivalent to 
adding vacancies respectively^!^. While in the latter, 
the lattice results from a spontaneous translational sym- 
metry breaking driven by the pressure as shown in the 
Fig. 2, if there are vacancies or interstitials in the ground 
state has to be self-determined by ground state energy 
minimization. There is usually no particle-hole symme- 
try for excitations. Due to this absence of symmetry, 
the number of vacancies is usually not equal to that of 
interstitials. So the theory developed in this paper on 
bulk *He is different from the lattice theory developed 
jj,43,50 giiQ-y^rn the appendix A and i n^°i^^ , one 
common fact of both supersolids is that both are due to 
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vacancies or interstitials. Lattice supersolids can also be 
described by doping the adjacent CDW either by vacan- 
cies or interstitials, so are classified as two types SS — v 
and SS — i. In the hard-core limit, SS ~ v and SS — i 
arc simply related by P-H transformation. However, as 
shown in section IV, in IIelium-4 supersolid, there is no 
particle-hole symmetry relating Tss~v to Tss-i ( Fig. 3 ) 
! 5*0 the coupling constant g in Helium 4 plays a similar 
role as the chemical potential n in the lattice models, the 
gap A(p) in the NS-PH which tunes the distance from 
the NS-PH to the SF plays a similar role as the gap in 
the CDW which tunes the distance from the CDW to the 
superfiuid . It was shown in^, the lattice supersolids ex- 
isting at commensurate 1/2 filling factors in frustrated 
lattices such as triangular lattice is just the coexistence 



of SS-v and SS-i. Combined with the results in^, we 
conclude that '^He supersolid can exist both in bulk and 
on substrate, while although H2 supersolid may not exist 
in the bulk, but it may exist on wisely chosen substrates. 
Lattice supersolid could also be realized in optical lat- 
tices in ultra-cold atomic experiments^. However, in 
both continuum and on lattices, SS states could be un- 
stable against phase separations. For example, the va- 
cancies or interstitials in the in-commensurate solid can 
simply move to the boundary of the sample instead of 
boson condensation, namely, it will turn into a commen- 
surate solid. This case is included in the C-NS case in 
the paper anyway. Due to its negative compressibility, 
the instability of lattice SS against phase separation was 
demonstrated in some lattice models in^. 
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